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I. INTRODUCTION 



A guided missile is one which is usually fired in a direction 
approximately towards the target and subsequently receives steering 
commands from the guidance system to improve its accuracy. 

In recent years, the application of Bank-to-Tum (BTT) guidance to 
tactical missiles has generated conciderable interest, motivated by certain 
advantages that a Bank-to-Tum control configuration can offer. 

First, against high performance threats, there is a need for defensive 
missiles to develops increasingly higher lift accelerations, which requires 
banking maneuver to properly direct the control vector. 

Second, the need for greatly increased stand off ranges has led the 
design of tactical missiles towards airbreathing propulsion systems such as 
the ramjet. These designs generally have configuration geometries that are 
not cruciform, due to exposed inlets beneath the vehicle. As a result, there 
are often stringent limits on the sideslip angle that can be developed during 
engine operation, and this too, dictates some types of Bank-to-Tum control 
scheme. 

This thesis investigates the controllability of a Bank-to-Tum missile, 
when one or more control surfaces are damaged, at fixed flight conditions 
(i.e. constant altitude, Mach number and missile weight and inertias). 

The approach is to design a control reconfiguration module to be placed 
between the flight control laws and the control surfaces. In the face of 
missile damage the reconfiguration module would reallocate contml 
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commands to the unimpaired control surfaces so that flight critical pitch, 
yaw, and rolling moments would be preserved to the maximum extent 
possible. Using the reconfiguration module control laws, the existing flight 
control laws would not be altered. 

In the first part of the thesis the control surfaces conventions, 
definitions and the manoeuvre by banking are reviewed. 

The second and third part reviews the design of the longitudinal and 
lateral channels and Involves the uncoupled pitch, roll and yaw channels for 
the circular and elliptical airframes of a B.T.T. missile. The resultant 
continuous open loop designs are analyzed and found to be In accordance 
with the desired specified requirements. Applying the convention control 
theory, the corresponding modified longitudinal and lateral autopilots are 
obtained, analyzed and compared with the original system. The effect of 
damage to one or more control surfaces, on the controllability of each 
system, is obtained, analyzed and compared with the unimpaired models 
respectively. 

Finally, the combined linear open loop pitch and roll channel autopilots, 
for the circular airframe configuration, are Investigated. The effect of 
damage to one or more actuators, on the control of the combined unimpaired 
system Is obtained and analyzed. The most severe type of Impairment 
considered In this thesis, I.e., the hard over type. Is extensively Investigated 
In Chapter V. Corrective actions are then developed and presented. 

The robustness. of the complete system Is analyzed In chapter five. 
Conclusions and recommendations for future study are then given. 
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The Investigation and analysis In all the above cases was performed 
using the existing at Naval Postgraduate School (NPS), Optimal Systems 
Control Program (OPTSYS) for the continuous system. Robustness Design 
Program (POPLAR) for the robustness of the system, and EASYPLOT 
program for the plotting of different kind of Impairments and comparisons. 
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II. GUIDED niSSILE CONTROL SURFACES CDHVEHT IDNS AND HAMCEUVRE 

BY BANKING 



In this chapter, before going into mathematical details concerning the 
motion of a guided missile in space as a result of guidance commands, some 
definitions and discussion concerning the control surfaces are given. 

A. GENERAL 

The task of a Missile Control system is to manoeuvere the missile 
quickly and efficiently as a result of the following two control singnals: 

(a) Detection whether the missile is flying too high or too low, or too 
much to the left or right. 

(b) Measure the deviations or errors and develope commands to the 

control system to reduce these errors to zero. 

B. CONTROL SURFACES DEFINITIONS 

Suppose that the control surfaces 1 and 3, as is shown in the Figure 
2.1, are mechanichally linked together such that a servo must impart the 
same rotation to both. Then these surfaces are essentially elevators . 
Positive elevator deflection produces a negative force in the 
z-di recti on and an anti -clockwise moment about the y-axis. 

If now surfaces 1 and 3 each have their own servo it is possible for 
them to act as ailerons . Positive aileron deflection produces an anti- 
clockwise moment about the x-axis. 
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If looking In the direction y one surface Is rotated 6 ^ clockwise and 
the other surface 6° antl-clockwise then a pure couple Is imparted to the 
missile about the fore and aft axis and this will tend to make the missile 
roll. Such control surfaces are now called ailerons, and by doing the same 
thing to control surfaces 2 and 4, the turning effect of the ailerons can be 
doubled. 



->y 

1 

2 

Figure 2.1 Control surfaces looking from rear of the Missile. 




If the aerodynamics are linear I.e., the normal forces are proportional 
to incidence then the principle of superposition applies. Commmands for 
elevator rudder and aileron movements can be added electrically resulting in 
unequal movements to opposite control surfaces. Positive rudder 
deflection produces a positive force in the y-direction and negative 
moment about the z-axis. In this way we have the means to control roll 
motion as well as the up-down i.e., pitch motion and left-right i.e., yaw 
motion. 
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C. CONTROL SURFACE CONVETION 

As is shown in the Figure 2.1 deflections 6^, 62, 63, 64 ore defined os 
positive if the leoding edge of the control surfoce is rotated clockwise as 
one looks outwards along the individual hinge axis. The total aileron 
deflection is defined as: 

Aileron deflection 6 = 1/4 (61+62+63+64) (II.C- 1 ) 

or 1/2(61+63) or 1/2(62+64) (II.C- 2 ) 

if only two surfaces act differentially. 



D. 0ANK-TO-TURN CONTROL SURFACES CONFIGURATION 
The same kind of analysis can be applied for the Bank-to-Tum missile 
configuration, shown in Figure 2.2. Note that the arrow on the control 
surface Indicates the direction of positive leading edge deflection. 







Tail 


Tail Incidence 


y 


Command 


^1 


5 2 


6 3 


64 




* ®Pc 




+ 


- 




\ 




* ^Yc 


+ 




+ 




* ®RC 


+ 


- 


- 


+ 



Figure 2.2 B.T.T Aerodynamic Sign Convention. 

The maneuver plane is the iTg-ZB plane. The downwards movement is 
defined as positive. This is the configuration investigated in this thesis. 
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The B.T.T. missile configuration Is tall -control led using four Identical 
control surfaces. The control surfaces are located flush with the body base 
with a i30 degrees dihedral. With tall control configurations the pitching 
moment or control effectiveness Is quite high In view of the fact that the 
tall moment arm is generally fairly large. Also there Is no adverse 
downwash from deflection of the forward surfaces, in comparison with the 
canard control configuration. A more complete model of the B.T.T. circular 
confinguratlon Is also shown In Figures 3.1 and 3.4 , taken from Arrow, A. 
IRef. 11. 

E. MANEUVER BY BANKING. 

The Conventional method of altering course to the right or left In a 
missile or aircraft Is to use the ailerons to bank I.e, roll by an angle (p, as Is 
shown In Figure 2.3. If the lift force Is Increased slightly by the use of the 
elevators so that the vertical component of lift equals the weight then 
there Is a horizontal component of lift equal to the total lift times sin 0. It 
Is this component of lift which causes the flight path to change. An exact 
analysis of the maneuver Is not simple as the airflow over the outer control 
surface Is faster than that over the Inner control surface due to the rate of 
turn. In practice a small amount of rudder Is applied In an attempt to make 
the general airflow directly along the fore and aft axis of the cordinated 
B.T.T. missile and In the plane of the control surfaces. In this condition 
there Is no “sideslip* and hence no net sideforce. This Is the preferred 
method of manoeuvring since lifting forces are most efficiently generated 
perpendicular to the control surfaces; the llft-to-drag ratio Is a maximum 
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In this condition. Even so there is some additional induced drag when 
manoeuvring. An additional manoeuvring force Is obtained by Increasing the 
Incidence and this means that the normal force Is inclined backwards to the 
velocity vector by the angle of incidence. There Is a useful component equal 
to the normal force times the cosine of the angle of incidence and a 
component proportional to the sine of the incidence opposing the velocity 
vector. This "Induced drag" will always occur whenever a lifting force is 
produced irrespective of the method of manoeuvre. 
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Figure 2.3 Forces Acting on a Banking Missile and Coordinated Missile 
Motion for coordinated B.T.T. Control Policies. 



Because it is desired that the B.T.T. missile avoid negative angles of 
attack, the missile is forced to roll about its velocity vector when the 
desired manouvre direction is in the negative angle-of-attack direction. 
Although some negative angles of attack may be possible with certain 
missiles. 
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Coordinated motion or zero sideslip angle Is achieved by directing the 
body fixed pitch axis of the missile at the missile velocity vector so that 
there is no component of missile velocity vector along the body fixed yaw 
axis of the missile. The Figure 2.3 (b) shows the attitude of the missile 
body with respect to its velocity vector vT When an upward command Is 
applied, I.e., (p = 0, the missile body moves upward with Its pitch axis 
directed at the velocity vector until it reaches manoeuver level or angle of 
attack. No roll motion Is required to maintain coordination of this 
manoeuver. The circular airframe configuration of this thesis Is a planar 
configuration which is symmetric about its wing planes and therefore have 
two perfect directions which are normal to the planes of Its wing. 
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III. LONGITUDINAL MODEL hETHODQLOGV 



A. INTRODUCTION 

In this Chapter the design of the continuous uncoupled pitch channel 
(longitudinal) autopilot Is reviewed, as was analyzed by Arrow, A. and 
Karalskos, I. [Ref. 1 and- 2). A modified longitudinal autopilot Is then 
designed with the four control surfaces split Into Independent components, 
for the circular configuration model shown In Figure 3.1. 

Then a comparison between the two systems (original and modified) 
shows that are Identical. The transient response plots and their analysis 
demonstrate the performance of the modified pitch channel. 

Having established the validity of the modified system one then uses 
the multiple-input multiple-output control theory and state-space 
representation [Ref. 3], to Investigate the performance of the modified pitch 
channel autopilot for various classes of damages. 

For this thesis, two classes of Impairments are considered namely, 
locked or hard over and Inoperative with the control surfaces centered. The 
plots of a modified pitch autopilot without Impairment are used as a basis 
for comparison, the Impairment Induced effects on normal acceleration, rate 
and tall Incidence are assessed. 

A general block diagram, of a B.T.T. autopilot Is shown In Figure 3.2 
taken from [Ref. U Inertial acceleration commands are applied In polar 
coordinates (I.e., magnitude of the command ( % ) applied to the pitch 
autopilot and the direction ( 9c ) is applied to the roll autopilot. Achieved 
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manoeuver plane or Inertial acceleration In rectangular coordinates ( I.e., i\i 
and tty) Is determined by resolving achieved body-fixed accelerations ( I.e., 
ttzand tty) through missile roll angle (9). 

Before going Into the details concerning the various kinds of damages 
It Is prefered to make here clear some terms used In this work. 

ORIGINAL SYSTEM is called the unmodified pitch channel 
autopilot presented in [Ref. 1]. 

UNIMPAIRED SYSTEM Is called the modified, I.e., with the split 
control surfaces system, and without any kind of damage. 

IMPAIRMENT, HARD OVER OR INOPERATIVE are terms expressing 
a specific kind of damage to one or more control surfaces. 



6. UNCOUPLED LINEAR PITCH CHANNEL (ORIGINAL) 

1. Aerodunamlc Model 

The linear pitch channel dynamic model which Is used for the 
circular airframe. Is shown In Figure 3.3, with the following requirements : 

a. In actuator command branch, the high attenuation frequency must be > 
15 db at 100 rad/sec and zero angle-of-attack. This requirement 
limits the autopilot speeds of time response. 

b. Relative stability with gain margins > 6 db, phase margins > 30 

degrees with a goal at 1 2 db and 50 degrees. 

c. An acceleration time responce with 6Z% time constant of 0.5 seconds 
for a step command of acceleration at the flight condition of Interest, 
overshoot < 10%, and zero steady-state error to reduce variations of 
guidance navigation gain. 
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The linear pilch channel dynamic model which Is used. Is shown In 
Figure 3.3. This Figure is token from [Ref. 11. The model is linearized for 
stability studies and the following three assumptions were made : 

a. Manoeuver plane is the plane ^ , as is shown in Figure 3.4 In 
which the arrow Indicates the direction of positive leading edge 
deflection. 

b. The missile is trimmed in pitch, at fixed values of angle of attack (a), 
pitch angular rate (q), and pitch tall incidence (5p), that is My= 0. 

c. The roll rote (P) of the missile is constant. 

Using the foregoing assumptions, and the approximation that the 
variation in forward velocity (u), is equal to zero a simplified set of 
longitudinal equations of motion are constructed applicable to the short 
period of motion. These equations describe the two degree of freedom, short 
period longitudinal mode. 

Linearized aerodynamic derivatives for M = 3.95 and Altitude = 60 
KFT are provided in Table 3.1, as is given in (Ref.ll. 

2. Pitch Control Law 

The longitudinal control low model for the circular airframe was 
developed from the block diagrams and information provided in (Ref. 1], and 
is shown in Figure 3.5. Lead-log filters were used in the system to shape 
the output response, provide adequate gain and phase margins, and finally to 
prevent guidance noise saturation problems. The guidance requirement of 
zero steady-state error is satisfied by using on integrator in the 
acceleration error branch of the circulare control low. The role error 
conpensotion determines the high frequency attenuation and was used to 
minimize the effect of aerodynamic variations on acceleration time 
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response. The acceleration error compensation determines the acceleration 
time response. 

The blocks In the pitch control law of Figure 3.5 for the circular 
model represent the basic transfer functions In which a normal acceleration 
command (hze) equal to 1 gee. Is applied to the pitch control low which uses 

measurements of: 

a. Missile body pitch angular rate (q) and 

b. Pitch normal acceleration (n^) 

to determine the required actuator comand (6pc>. The actuator Is modeled as 

a first order lag at 188.4 rad/sec and Is shown In Figure 3.2. 

3. Basic Transfer functions of the Aerodynamic Model 

The equations which are represented by the block diagram of the 
aerodynamic model, of the uncoupled pitch channel are presented below and 
taken from lRef.21 . 

a. The pitch aerodynamic transfer function for the pitch 
angular rate (q) about the Yg axis: 

q -1.3361s- 0.15 

— = ( 1 /sec) (1 1 1.8.3- 1 ) 

6p 22.545 X 10-3 s2 ♦ 3.3634 X 10-3 s + 1 

b. The pitch aerodynamic transfer function for the achieved 
maneuver normal acceleration (t^) In the ^ direction: 

n^ 106.074 X 10-3 s2- 18.82 

= (gs/rad) (III.B.3-2) 

6p 22.545 X 10-3 s2 ♦ 3.3634 X 10-3 s ♦ 1 
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Figure 3.1 B.T.T. Model of Circular Configurolion; Scale 1/6. 
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Figure 3.4 Aerodynamic Sign Convention and Axes. 
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TABLE 3.1 



LINEARIZED AERODYNAMIC DERIVATIVES (M=3.95 ) 
FOR CIRCULAR AIRFRAME 





a = 0 <> 


a = 10 o 


a = 20 ® 




- 0.065 


- 0.082 


- 0.1 1 1 




- 0.025 


- 0.019 


- 0.003 




0 


- 0.009 


- 0.020 


Cy 6 Y 


0.021 


0.022 


0.028 


Cnfi Y 


- 0.050 


- 0.053 


- 0.062 


ClsY 


0 


- 0.016 


- 0.038 


Cysr 


0 


- 0.009 


- 0.022 


CnfiR 


0 


0.018 


0.044 


^ISR 


0.031 


0.035 


0.044 




0.15 


0.17 


0.22 




0.04 


0.04 


0.05 




- 0.060 


- 0.065 


- 0.118 




- 0.080 


- 0.095 


- 0.115 



Reference Center of Gravity (C.6) at 0.6 1 
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4. Analysis for Circular Airframe 

The analysis of the uncoupled pitch channel autopilot is based on 
the time responses and Bode plots of maneuver plane acceleration, body 
angular rates and tall Incidence angles, obtained by 0PTSY5 control program 
of Naval Postgraduate School (NPS), which requires the differential 
equations to be In the state-space variable format: 

X = F X ♦ G U (III.B.4-1) to 

define the system. In this system the F matrix is the open loop dynamics 
matrix (system or plant) and the 6 matrix is the control matrix. The 
variable assignements are X as the state vector and U as the control Input 
vector. It follows that X Is the time derivative of X. 

Using state-space representation of a system. In which the 
forcing function Involves derivatives, the following equations are derived 
and can be modeled as a tenth order system, for the uncoupled pitch 
channel : 



X,= -150 X, +150 Xfi + 705.75 6 p 


(III.B.4-4) 


Xj = X 3 - 1.3531 nzc ♦ 1.3531 X, 


(III.B.4-5) 


X 3 = -5 X 3 + 6.572 nzc - 6.572 X, 


(III.B.4-6) 


X 4 = X 5 - 59.2637 8 p 


(III.B.4-7) 


Xg = -H9.186 *10-3 Xg -44.3557 X^ + 1.7887 8 p 


(III.B.4-8) 


Xg = X 7 - 701.92x10-3 5p 


(III.B.4-9) 


X 7 = -149.186x10-3 X 7 - 44.3557 Xg - 1043.36 5p 


(II 1.8.4- 10) 


8 p = -188.4 5p + 188.4 8 pc 


(III.B.4-1 1) 



39 



= -0.H3 5po -2.2308 Xj + 2.2308 X 4 - 340.58x10-3 X 3 
♦ 460.839x10-3 - 460.839x10-3 X, + 340.58 

X 1 0-3 Xg - 20. 1 84 5p ( 1 1 1.8.4- 1 2) 

riz = -( X 7 - 887.124 6 p + 886.422 6 po ) (1 1 1.8.4- 13) 



where the state variables are: 

• Xi X 2 : pitch angular rate state variables 

• Xg X 7 : achived body acceleration state variables 

• Xg . output of acceleration filter 

• Spc : input command in actuator network 

• Sp . tail incidence angle 

Expressing the above equations as a set of first order matrix 
in the format of (lII.B.4-1) the plant and control matrices are: 



F= 



-150 


0 


0 


0 


0 


0 


0 


0 


0 150 


1.353 


0 


1 


0 


0 


0 


0 


0 


0 


0 


-6.572 


0 


-5 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 ■ 


-59.26 


0 


0 


0 


0 


0 


-44.36 


-0.15 


0 


0 


1.789 


0 


0 


0 


0 


0 


0 


0 


0 


1 - 


0.7019 


0 


0 


0 


0 


0 


0 


0 - 


44.36 


-0.1492 


-1043 


0 


0 


0 


0 


0 


0 


0 


0 


0 


-188.4 


188.4 0 


-0.461- 


■2.231 


-0.3406 2.231 0.3406 


0 


0 


-20.18 


-0.143 0 


0 


0 


0 


0 


0 


0 


-1 


887.1 


-886.4 0 



(III.B.4-14) 



The control matrix 6 is following: 
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6 = 



0 

-1.3531 

6.572 

0 

0 

0 

0 

0 

0.4608 

0 



(III.B.4-15) 



Where: U = control input to the system. (III.B.4-16) 

According to the [Ref. 1], using a step input function which 
represents the ‘1 gee Command’ and a trim angle of attack a«=0, the time 

response plots of the state-variables of interest were obtained from the 
OPTSYS program and are represented at the Figures 3.6, 3.7 and 3.8. Typical 
performance criteria used to characterize the transient response to the one 
gee unit step input include overshoot, delay time; rise time, and settling 
time. The above-mentioned criteria are defined with respect to the step 
response; 



• tidximum Overshot The maximum overshoot is defined as the largest 
deviation of the output over the step input during the transient state. 
The amount of maximum overshoot is also used as a measure of the' 
relative stability of the system. The maximum overshoot is 
represented in this thesis as a percentage of the final value of the 
step response: that is. 



max. overshoot 

percent maximum overshoot = * 1 00% (III.B.4- 1 7) 

final value 



41 



• Delay time. The delay time Is defined as the time required for the step 
response to reach 50% of Its final value. 

• Rise time. The rise time is defined as the time required for the step 
response to rise from 10% to 90% of its final value. Frequently the 
rise time is also refered to as the time constant of the system. 

• Setiiing The settling time is defined as the time required for the 
step response to decrease and stay within a specified percentage of 
Its final value. For this thesis Is specified arbitrarily at 0.01%. 

The four quantities defined above give a direct measure of the 

transient characteristics of the system to the step Input: 

a. The pitch normal acceleration (n^ ro N^), Figure 3.6, g's vs. 

time which has: 

( 1 ) A 0.5 rise time or time constant. 

(2) A minor 0.005 steady state error and 2.2% overshoot. 

(3) A 0.28 sec delay time and 2 sec settling time. 

b. The pitch angular rate (q or Q), Figure 3.7, rad/sec vs. time, 

which has: 

(1) The maximum body angular rate (*0.17 rad/sec) occurs at 0.28 second 
when the tall Incidence Is at 0.02 rad angle of deflection. 

(2) From the moment at 0.28 sec the rate increases and reaches at 2 sec 
its steady-state value (-0.01 rad/sec). 

(3) The -0.01 rad/sec steady-state constant rate Is due to the fact that 
one has a type zero system and therefore has a finite error to a step 
input. 

c. The pitch tail incidence (Sp or DP), Figure 3.8, rad vs. time, 

which has: 

(1) Actuator deflection angle 0.05 rad. 
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(2) Initially and within the 0.28 sec delay time a small transient occurs 
with an insignificant fluctuation. 

(3) An overshoot \.Z% and a 2.0 sec settling time. 

The above analysis has revealed that the perfomance of the system is 
according to the previous mentioned requirements. 

Also should be mentioned here that according to the Eigenvalues found 
executing the OPTSYS program, the continuous open loop system is stable, 
since the s-plane poles are: 



s, = -159.724 


♦ j 18.9921 


S 2 = -159.724 


- j 18.9921 


S 3 = -8.29048 


♦ J8.05932 


S 4 = -8.29048 


- J8.05932 


S 5 = -3.75925 


♦ j2.51463 


Sg = -3.75925 


- j 2.5 1463 


S 7 = -0.14393 


♦ jO 


Sg = -0.00000425 + jO 


Sg = -0.04593 


+ J6.65959 


S|Q- -0.04593 


- J6.65959 



5. Design of the Modified Uncoupled Pitch Autopilot 

To study reconfigurable flight controls the missile model should 
be capable of using the set of four available control surfaces to produce the 
required forces and moments. 

To achieve this each control surface ( S^, 62 , 63, S4 ) is considered 
as an independent component. 
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The equations are then coupled so that a component of control 
surfaces used as a single or pair, will produce the appropriate moments. For 
example 6t deflected at an certain angle, while the other three remained at 
a zero angle, will produce primarily a rolling moment, and deflected 
together with 82 no roll effect is present, as we will see in Chapter V. 

Taking under consideration the convention control surfaces theory 
of Chapter II., for the e.T.T. missiles, and for the configuration shown in the 

Figure 3.4, one can easily show, for the continuous pitch autopilot that the 
input command in the actuator network (fip^) con be split up into independent 

components, using the following equation which compute the deflection 
inputs to Equation ( III.B.4-12 ). (Also refer to Figures 2.2 and 3.4 which 
show the control surface positions and corresponding positive deflections). 

*6pc = (-61 +82 -83 +84 ) / 2 (III.B.5- 1 ) 

Replacing the input vector in Equations (III.B.4- 1 1 ), (III.B.4-12) 
and (III.B.4- 13) with the above Equation and applying the superposition 
principle, gives the following modified airframe equations: 

1 

8p = -188.48p ♦ 188.4 — (- 8, + 82 - 83 ♦ 84) (III.B.5-2) 

2 

1 1 

— (- ♦ 82 - 63 ♦ 64) = -0.143 — (-8t ♦ 82 - 83 + 84) - 2.23O8X2 

2 2 

♦ 2.23O8X4 - 340.58x10-3 X3 + 460.839x1 0-3n2c 

- 460.839x1 0-3X| ♦ 340.58x10-3X5 - 20.1848p 

(III.B.5-3) 
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1 

rizo = - (X 7 - 887.124)5p + 886.422 — (-5, + 62 - 63 + 64 ) (III.B.5-4) 

2 

In Figure 3.6 the modification of the commanded tail incidence is 
presented. 




(a) 




Figure 3.6 (a) Original System; (b) Modified System. 



Rearranging the above equations and using the transfer functions 
of the system for the uncoupled pitch channel. Equations (III.B.4-4 to 
III.B.4-13), the following equations are derived and can be modeled in a 
thirteenth order system for the same channel: 



X, = -150Xi ♦ ISOnz 

X 2 = X 3 - 1.3531nzc+ 1 . 3531 X 1 

X 3 = - 5 X 3 + 6.572nzc - 6.572X| 



(III.B.5-5) 

(III.B.5-6) 

(III.B.5-7) 



X 4 = Xs - 59.26376p (III.B.5-8) 

X 5 = - 149. 1 86X 1 0 - 3 X 5 - 44 . 3557 X 4 + 1 .78875p (III.B.5-9) 

Xg = X 7 -70 1 .92X 1 0-36p (III.B.5- 1 0) 

X 7 = - 1 49. 1 86 x 1 0 - 3 X 7 - 44.3557Xs- 1 043.366p (III.B.5- 1 1 ) 

^ = -188.46p - 94.26t ♦ 94.262 “ 94.263 + 94.264 (III.B.5- 12) 



= -0.1436| ♦ 1.1 154 X 2 - 1.1 154 X 4 + 0 . 17029 X 3 
-0.23042nzc + 0.23042Xi - 0.17029Xs + 10.0296p (III.B.5- 13) 

^2 = -0.14362 - 1.1 154X2 ♦ 1*^ 154X4 - 0.17029X3 

♦ 0.23042nzc - 0.23042X2 + 0.17029Xs - 10.0926p (III.B.5- 14) 

^3 = -0.1436| + 1.1 154 X 2 - 1.1 154 X 4 ♦ 0 . 17029 X 3 

- 0.23042nzc ♦ 0.23042X, - 0.17029Xs + 10.0296p (III.B.5- 15) 

^4 = -0.14364 - 1-1 154X2 + 1.1 154X4 - 0.17029X3 

♦ 0.23042nzc - 0.23042X| ♦ 0.17029Xs - 10.0296p (III.B.5- 16) 

fiz = -X 7 ♦ 887.1246p ♦ 443.21 16| -443.21 162 + 443.21 163 

- 443.211 64 (III.B.5- 17) 

We thus have a multIple-input-muUlple-output (hIMO) system 
with the following order 

• X : 13x1, the time derivative of X matrix. 

• F : 13x13, the open loop dynamics matrix or plant. 

• X : 13x1, the variables assignements matrix. 

• 6 : 13x5, the control distnbutlon mathx. 

• U : 5x1, the control Input matrix. 

The above plant control and output matrices are presented at the 
Appendix A of this thesis. 
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6. Anali^sis and Comparison with the Original System 

Applying the same step Input function as In section III.B.4, which 
represents the "1 gee command', and zero degrees deflection command for 
the four control surfaces (i.e., ^20 ^30 ^4c = response plots 

of the following state-variables of Interest are obtained and presented In 
the Figures 3.10 through 3.13. 

The performance of the system ( q, 8p ) os shown in Figures 

3.10 to 3.12 Is the same as that shown In Figures 3.6 through 3.6. 

In Figure 3.13 Is shown the time response of the control surfaces 
deflection angle, for the Input of an acceleration equal to 1 gee command. 

Deflections ( 82 * ^31 ^4 or D"s ), Figure 3.13, rad vs. time, 

which has: 

(1) The 8 t and 83 have the same sign (negative), as It was expected 
according to the sign convention of Figure 2.2. 

(2) The 82 and 84 , have the same response as above but with a positive 
sign. 

(3) Insignificant overshoot. 

(4) A 2 seconds settling time. 

In the same Figure 3.13 the deflection of the control surfaces Is 
Illustrated. The presented deflection Is the required position of each 
surface In order to achlve the commanded one gee acceleration (downwards). 

We procceed now to the analysis of the Impaired control surfaces 
for the modified system. 
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PITCH NOHMAL ACCELEUATION,NZ 

0.0 0 I 0 2 0.3 0.4 0 5 0.6 0.7 0.8 




0.0 0.6 1.0 1.5 Z.a Z.l 3.0 



TIME - SEC 



Figure 3.7 Pitch Normal Acceleration vs. Time; Uncoupled Pitch Channel; 
Continuous Open Loop System; a« = 0; Step Input = 1 g's. 
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Q, PITCH ANGULAR RATE-RAD/SEC 

- 0.175 - 0.150 - 0.125 - 0.100 - 0.075 - 0.050 - 0.025 0.000 
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C. CONTROL EFFECTIVENESS DUE TO THE CONTROL SURFACES 

IMPAIRMENT AND COMPARISON 

The time response plots of the system are obtained In this section 
for two types of damages. The first type Is the locked or hard over 
surface at a small angle and the second type Includes the Inoperativlty of 
one, two or three actuators. The second type of damage could also 
represent the destruction of the control surface(s), I.e., shot away. 

Note that the first' phase of the design approach for the CBTT 
autopilot was to design each channel Independently with all coupling 
between channels removed. Sufficient high frequency attennuatlon was 
added for actuator and missile design so that the resulting missile body 
angular rates and control surfaces motion would represent a practical 
missile design. A relationship was established among the relative speeds 
of response of the uncoupled channels In order to meet CBTT requirements. 
The cross-coupling due to control surfaces Impairment between the 
longitudinal and lateral motion of the missile is ignored when each control 
channel Is treated Independently. Cross-coupling of the controls Is thus 
considered In Chapter V. 

1. Hard Over One Control Surface 

One control surface (6^), was selected to be deflected or hard 
over, 1.e., It remained at a fixed angle while the other three control 
surfaces have a variable deflection angle. 

The control Input was a step Input of one gee In the normal 
acceleration as In the original and unimpaired models. For the damaged 
control surface (6,) a hard over deflection of 0.2 rad was considered. The 
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time response plots of the state-variables of Interest were obtained and 
analyzed with the performance criteria as defined In section (III.B.4): 

a. Pitch normal acceleration (n^ or N^), Figure 3.14, g's vs. 

time, which has: 

(1) Initially and with a fast response, the missile body tends to go 
upwards reaching a 0.6 gee acceleration. This Is caused by the fact 
that at the beginning of the command the damaged surface was 
deflected at a *0.Z rad angle of deflection while the other three 
were at their steady position. Then the missile goes downwards 
until the -1.25 gee acceleration is achieved. This is due to the 
deflection of the unimpaired surfaces which are commaned In such a 
way as to counteract the damaged (St) surface. From this point (0.2 
sec) the acceleration rises constantly In order to achieve the one 
gee command. 

(2) A rise time 0.75 seconds and a negligible overshoot. 

(3) A delay time about 0.5 sec 

(4) Steady-state error 0.005% which Is Insignificant. 

b. The pitch angular rate (q or Q), Figure 3.15, rad/sec vs. time, 
which has: 

(1) Initially and due to the fact that St Is locked at a positive angle of 
0.2 rad, to achieve the acceleration response, the required body 
angular rate the first tenth of the second Increases very fast to 
0.48 rad/sec. Then It begins to decrease crossing the zero rate at 
0.2 sec, corresponding to the max decceleratlon (-1.25 gee's). Then 
It changes direction for the next 0.2 sec until the value of -0.36 
rad/sec Is achieved. From this point the rate Is increased again 
until the constant value of -0.01 rad/sec Is achieved at the 
settling time 1.4 sec. 

(2) The Initial Increase of the rate Is due to the produced pitch moment 
caused by the deflection of the damaged control surface while the 
other three are at their steady position. 
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c. Pttch tall Incidence (6p or DP), Figure 3.16, rad vs. time 



which has: 

(1) To achieve the acceleration response and due to the fact that is 
locked at an angle of 0.2 rad, initally more control surface 
deflection is required ( -0.13 rad). At 0.2 sec, when the 
acceleration begins to increase and the angular rate is zero, the tail 
incidence has the positive angle of 0.03 rad and a fluctuation begins 
before being stabilized at a constant angle of 0.06 rad. 

(2) Settling time equal 'to 1.8 sec. 

d. Deflections ( S^, S2, S3, 54 or D's ), Figure 3.17, rad vs. time. 

This Figure illustrates the time response of the deflection 
of the four control surfaces,when one of them (S|) is locked at an angle of 
0.2 rad, in order to achieve the acceleration response. Initially while S2, 
S3 and S4 control surfaces have zero angle of deflection, S^ is already at 
the 0.2 angle. Then S2 and S4 are deflected at a positive angle 0.17 rad and 
S3 at a negative 0.17 rad, according to the sign convention shown in Figure 
2 . 2 . 



A comparison between Figure 3.17 and 3.13 of the unimpaired 
system, shows that to achieve the acceleration response the impaired 
system requires more control surface deflection. 

Finally to cover completely this class of impairment, the case of 
a negative angle (-0.2 rad) at the S^ surface is examined. The time 
response of the same state variables are obtained and analyzed. 

a. Pitch normal acceleration (n* or Nz), Figure 3.18, g's vs. 

time, which has: 

(1) Initially and with a fast response, the missile body tends to go 
downwards reaching a 0.6 gee acceleration. This is caused by the 
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fact that at the beginning of the command the damaged surface was 
deflected at a -0.2 rad angle of deflection while the other three 
were at their steady position, then the missile goes upwards until 
the +1.75 gee acceleration Is achieved. The produced overshoot Is 
due to the deflection of the unimpaired surfaces which are 
commanded In such a way as to counteract the damaged (6|) surface. 
From this point (0.3 sec) the acceleration deacreases constantly In 
order to achieve the one gee command. 

(2) A small 0.18 second rise time (reason for the big overshoot). 

(3) A delay time about 0.1 sec 

(4) Steady-state error 0.005%. 

b. The pitch angular rate (q or Q), Figure 3.19, rad/sec vs. time, 
which has: 

(1) Initially and due to the fact that 6) Is locked at a negative angle of 
-0.2 rad the required body angular rate the first tenth of the second 
decreases very fast to 0.68 rad/sec. Then It begins to Increase 
crossing the zero rate at 0.3 sec, corresponding to the max 
decceleratlon (1.75 gee's), and It changes direction for the next 0.5 
sec until the value of 0.13 rad/sec Is achieved. From this point the 
rate Is decreased again until the constant value of -0.01 rad/sec Is 
achieved at 1.8 sec settling time . 

(2) The big overshoot of the rate Is due to the produced pitch moment 
caused by the deflection of the damaged control surface while the 
other three are at their steady position. 

c. Pitch tall Incidence (8p or DP), Figure 3.20, rad vs. time 
which has: 

(1) To achieve the acceleration response and due to the fact that Sy Is 
locked at a negative angle of -0.2 rad, Initally more control surface 
deflection Is required (about 0.16 rad). At 0.3 sec, when the 
acceleration begins to decrease and the angular rate Is zero, the tall 
incidence has the positive angle of 0.01 rad. A fluctuation begins 
before the deflection being stabilized at a constant angle of 0.06 
rad. 
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(2) Settling time equal to 1.8 sec. 

d. Deflections ( 8 ^ 62. 5-4 or D's ), Figure 3.21, rad vs. time. 

This Figure Illustrates the time response of the deflection 
of the four control surfaces, when one of them (6^) is locked at a negative 
angle of -0.2 rad. In order to achieve the acceleration response. Initially 
and when the S2, S3 and S4 control surfaces have zero angle of deflection 
St Is already at the -0.2 angle. Then S2 and S4 are deflected at a negative 
-0.09 rad angle and S3 at a positive 0.09 rad. 

A comparison between Figure 3.21 and 3.17 indicates that when 
the surface S^ Is hard over at a negative angle: 

(1) The undamaged control surfaces In order to counteract the failure 
and achieve the commanded acceleration, are deflected at a smaller 
angle. 

(2) The undamaged control surfaces are commanded to turn at angles 
with opposite sign comparing with those In Figure 3.17 (S^ 
positively hard over). 

2. Damage One Control Surface 

The St control surface was selected to simulate the Inoperative 
actuator. The corresponding state variable St was therefore set equal to 
zero. The step Input of one gee remained the same. The Initial condition on 
the control surface was set at zero. This type of damage could also 
represent the distructlon of the control surface. 

For the above case the time response plots of the state variables 
of Interest were obtained and analyzed according to the defined 
performanse criteria In Section (III.B.4). Later, In Section 5, a comparison 
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DP, PITCH TAIL INCIDENCE-RAD 

- 0.150 - 0.125 - 0.100 - 0.075 - 0.050 - 0.025 0.000 0.025 0.050 0.075 
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Q, PITCH ANGULAR RATE-RAD/SEC 
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Figure 3.21 Control Surfaces Deflections vs. Time; Hard Over 6 i Control 

Surface; ag=0; Step Input; rijc=1 gee, 6ic=-0.2 rad, S 2 , 631 ^ 4 =*^- 
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between the different types of damages and the unimpaired system was 
performed. 

a. Pitch normal acceleration (n^ or Nz)» Figure 3.22, g's vs. 
time, which has: 

(1) Time constant or rise time 0.41 seconds. 

(2) A 0.005% steady state error. 

(3) Minor overshoot and delay time 0.35 sec and 

(4) A 1.6 sec settling time at which the acceleration response is 
achieved. 

b. Pitch angular rate (q or Q), Figure 3.23, rad/sec vs. time, 
which has: 

(DA maximum body angular rate of 0.17 rad/sec is achieved at 0.3 
second. 

(2) Minor overshoot and a settling time equal to 1.6 after which the 
angular rate remains constant at -0.01 rad/sec. This is due to the 
fact that the transfer function q/6p (Equation III.B.3-1) is a type 
zero system and has a finite error to a step input. 

c. Pitch tail incidence (5p or DP), Figure 3.24, rad vs. time. 



which has: 

(1) A 0.054 rad command from the actuator network. 

(2) A 0.4 sec delay time and negligible overshoot. 

d. Deflections (S^, h* ^4 o'" Figure 3.25, rad vs. time. 



which has: 

(DA positive deflection angle of 0.035 rad for the $2 and 64 actuators 
while the is deflected at a negative and equal angle 
(-0.035 rad), and 6^ as damaged remains constantly at zero. 
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(2) A 1.6 sec settling time for each of the unimpaired surfaces. 

(3) A rise time equal to 0.7 sec and minor overshoot correspondingly 
for the undmaged surfaces. 

(4) The deflection angle of the unimpaired control surfaces has a 
pitching moment effect on the system. 

A comparison between the deflection angle in Figures 3.13, 3.17, 
3.21, and 3.25 indicates that for the same step input: 

(1) When one control surface is Inoperative or shot away the system 
commands the undamaged surfaces at larger angles than those 
obtained with the unimpaired system. 

(2) When one control surface is inoperative or shot away the system 
commands the undamaged surfaces at smaller angles (in comparison 
with the hard over type of damage). 

3. Damage Two Control Surfaces 

In this case two control surfaces, on the same side of the 
missile were considered damaged and S 2 ). The two state variables 
and $2 therefore set equal to zero. The step input of one gee 

remained the same. The initial condition on the control surfaces was set 
at zero. Then the time response plots of the state variables were obtained 
and analyzed. 

a. Pitch normal acceleration. Figure 3.26, which has a 0.6 rise 
time, 0.38 sec delay time and 0.005% steady state error. Note the faulty 
behavior at the start of the response (0 to 0.15 sec). Because the system 
Is a nonminimum phase system (i.e., a zero lies in the right-half s plane) 
the transient response starts out in the opposite direction to the input 
i.e., there is an inherent delay in the transient response of the system. 
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b. Pitch angular rate. Figure 3.27, which has a maximum rate 
or -0.16 rad/sec at 0.3 sec. Then increases until the constant negative 
rate (-0.01 rad/sec) is achieved at 1.3 seconds. 

c. Pitch tail incidence. Figure 3.26, which has a 0.8 sec rise 
time and 0.3 sec delay time. The command from the actuator network is 
0.054 rad which is achieved at 1.6 sec settling time. 

d. Deflections (6f, 82, 83, 84 or O's), Figure 3.29, in which 
the two damaged control surfaces (81 and 82) remain always at zero 
(damaged) and the two others opposite deflected at the angle of 0.053 rad 
angle after a lapse time of 1.6 sec. Note that for this system with two 
impaired control surfaces more angle of deflection is required (In 
comparison with the system having one actuator shot away). 

4. Damage Three Control Surfaces 

In the last case three of the control surfaces 8|, 82 and 83 were 
damaged or inoperative. 

With a step input of one gee's for the commanded normal 
acceleration as in the unimpaired and other type damaged models, the 
plots of the state variables of Interest were obtained and analyzed. 

Figures 3.30 through 3.33 Illustrate the response of the output 
variables. Figure 3.30 shows that the commanded one gee acceleration can 
still be achieved in the absense of three control surfaces (8^, 82, 83). 
However, some oscilation is evident into the response. The loss of the 
surfaces has produced a more lightly damped system as evidenced in 
Figures 3.30 through 3.33. The oscillatory and convergent behavior of the 
pitch rate and the undamaged surface (84), shown in Figures 3.31 and 3.33, 
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Illustrate how the surface attempt to "fill In" for the absence of the 
inoperative control surfaces. Specifically the response of the pitch normal 
acceleration (Figure 3.30) has: 

(1) A 1.01 sec time constant. 

(2) A 0.38 sec delay time and 

(3) A 3.0 sec settling time at which the steady-state value of the 
commanded acceleration Is achieved. 

A comparison with the response of the corresponding deflection 
of the unimpaired system (Figure 3.13) and the other inoperative cases, 
Figures 3.21, 3.25 and 3.29, Indicates that as the number of the 
Inoperative control surfaces Increases, the deflection of the undamaged 
actuatoKs) also Increases and more time Is required to achieve the 
steady-state value. 

5. Comparison with the Unimpaired System 

In this section a performance comparison of each type of 
Impairment with the undamaged system Is performed and an analysis of 
their corresponding time response plots Is presented. 

In order to use the N.P.S. EASYPLOT program the state variables 
had to be input as data. This was accomplished by writing a computer 
program ( Appendix D). 

With the above data Figure 3.34 Is obtained. This corresponds to 
pitch normal acceleration (nx or NZ) for the unimpaired, hard over St* ^nd 

Inoperative or shoot away one or more control surfaces, configurations. 

The hard over case shows greater excursion In the normal 
acceleration with a peak value of -1.25 gee's at 0.2 sec being achieved. 
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NZ, PITCH NORMAL ACCELERATION-G'S 
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Q.PITCH ANGULAR RATE-RAD/SEC 
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NZ, PITCH NORMAL ACCELERATION-G'S 
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Figure 3,29 Control Surfaces Deflections vs. Time; 

Inoperative 6^ and 63 Control Surfaces; 3^=0; 

Step Input; nzc =1 gee, S,c = -0.2, 53c, Sjc, S4C = 0 rad. 
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The Inoperative cases are all quite similar except for the last 
case where only one control surface Is available. In this latter case an 
oscillatory behavior Is evident. The performance criteria to the one gee 
command are summarized In Table 3.2. 

TABLE 3.2 

PITCH NORMAL ACCELERATION PERFORMANCE 

System 

Unimpaired- Hard Over- Inoperative Control Surfaces 









One 


Two 


Three 


Time Constant(sec): 


0.50 ... 


,....0.41. 


.0.60.., 


0.60.... 


1.01 


Delay Tlme(sec): 




.....0.38.... 


0.40 


Overshoot(3>): 


2.2%..., 


...0.01%. 


0.01%.., 


...0.01%.... 


,...0.01% 



Figure 3.35 compares the unimpaired system with the three 
cases of control surfaces impairments. The curves Include the following 
four cases: 

(1) Original system 

(2) Inoperative one control surface 

(3) Inoperative two control surfaces 

(4) Inoperative three control surfaces 

This Figure shows that as the number of the damaged control 
surfaces Increases the delay time Increases also from 0.28 sec for the 
unimpaired system to 0.41 sec for the system with three Inoperative 
control surfaces. The Increase also of the time constant from 0.5 sec 
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(unimpaired system) to about 1 sec (three Inoperative actuators), 
Indicates that as the number of the damaged control surfaces Increases 
more time Is required to achieve the steady-state value of the normal 
acceleration. Also observe the initially faulty response of the nonminimum 
phase system to the step Input which decreases as the number of the 
inoperative control surfaces increases. 

Figure 3.36, shows that to achieve the acceleration response the 
system with the hard over 8 ^ type of damage, requires more angular rate 
(q or Q), compared with the unimpaired and the other three damaged 
system cases. The produced pitching moment Is evident. Initially and 
because the Sf control surface is deflected at the 0.2 rad angle while the 
other three are at the zero position, the first tenth of the second the rate 
Increases from zero to 0.45 rad/sec. The rate then begins to decrease 
very rapidlly to a rate of -0.38 rad/sec at 0.4 sec. 

Figure 3.37 compares the unimpaired system with three cases of 
control surfaces Impairments. The curves include the following four cases: 

(1) Original system 

(2) Inoperative one control surface 

(3) Inoperative two control surfaces 

(4) Inoperative three control surfaces 

This Figure shows that as the number of the damaged control 
surfaces Increases, the angular rate of change Is reduced, and more 
pronounced fluctuations occur, 1.e., the system becomes more oscillatory. 
These fluctuations are connected with the pitching moment Induced In the 
missile by the failed actuators. 
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Finally Figures 3.38 and 3.39 compare the unimpaired system 
with four cases of control surfaces Impairments. The curves Include the 
following five cases: 

(1) Original system 

(2) Hard over one control surface 

(3) Inoperative one control surface 

(4) Inoperative two control surfaces 

(5) Inoperative three control surfaces 

Figure 3.36 and 3.39 show that the required (steady-state) pitch 
tall Incidence for the uncoupled longitudinal autopilot is 0.054 rad 
(unimpaired or damaged). 
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NZ, PITCH NORMAL ACCELERATION-G’S 
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Figure 3.37 Pitch Angular Rate; Uncoupled open loop Channel: 
Unimpaired and Damaged Control Surfaces. 
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Figure 3.38 Comparison of the Pitch Tail Incidence between the 
Unimpaired System and the Four Cases of Impairment. 
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IV. LATERAL MODEL flETHODOLQGV 



A. INTRODUCTION 

In this chapter the design and analysis of the continuous uncoupled 
lateral (i.e., yaw and roll) channels of a B.T.T. autopilot derived in (Ref. 1 
and 4] is reviewed. A modified lateral autopilot is then designed with the 
four control surfaces split into independent components. The model of 
elliptical (yaw channel) and circular (roll channels) configurations are both 
shown in Figure 4.1. Both airframes have the same cross sectional area 
distribution and specifically, the circular cross sectional body has a closure 
ratio Abase/ Amax of 0.69 with Amax occuring at 68X missile body, whereas 
the elliptical airframe has a 3:1 cross section. 

A comparison between the two systems (original and modified) shows 
that they are identical. The transient response plots and their analysis 
demonstrate the performance of the original and modified yaw and roll 
channels. 

Having established the validity of the modified system one then uses 
the multiple-input/multiple-output control theory and state-space 

representation [Ref .31 to investigate the performance of the modified 
lateral autopilot for various classes of damages. 

For the investigation of the lateral channel, two classes of 
impairments at the control surfaces are considered: 

( 1 ) Locked or hard over and 

(2) Inoperative with the control surfaces centered or shot away. 
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The plots of the unimpaired modified yaw and roll autopilots are used 
as a basis for comparison. The Impairment induced effects on roll angle, 
rate and tall Incidence are assessed. 

B. UNCOUPLED LINEAR YAW CHANNEL (ORIGINAL) 

1. Aerodynamic Model. Airframe Configuration and Control Lew 

A general block diagram of a 6.T.T. autopilot with all Its channels 
Is shown In Figure 3.2. Inertial acceleration commands are applied In polar 
coordinates, i.e., magnitude of the command (ric) applied to pitch autopilot 
and the direction (^q) is applied to the roll autopilot. The yaw autopilot is 

slaved to the roll autopilot In order to minimize sideslip angle by 
coordinating the missile yaw and roll motion. Achieved maneuver plane 
accelerations in rectangular coordinates, (riz) and (r|y)^ were determined by 
resolving achieved body-fixed accelerations, (Hz) and (riy) through the 
missile roll angle (*p), i.e., the Euler angles 6 and ip are assumed to be 
sufficient small, so that the sideslip angle (B) Is minimized and the motion 
between the yaw and roll channels Is coordinated. 

A general diagram of the uncoupled yaw channel where the 
aerodynamic and yaw control law are Involved Is shown in Figure 4.2. In 

Figure 4.3 (Ref. 1) the yaw control law Is presented. The yaw normal 
acceleration (ir\yc) is not used to command the coupled 6.T.T. autopilot, 

although. It Is used for the design and analysis of the uncoupled channel. The 
command used by the coupled system Is shown In dashed lines and Is a yaw 
angular rate command (r^). The yaw control low is governed by missile body 
yaw angular rote (r) and yaw normal acceleration ( ity )• 
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The night conditions are the same as was mentioned for the pitch 
channel, (I.e., 60 Kft altitude and Mack number 3.95). The yaw control law 
determines the required yaw toil Incidence (6vc) command to the actuator 
network which Is approximated as a first-order lag at 30 Hz. 

2. B.equtrem.eD,tUii r .tri e .Cla.s.slMl Peslgg 

Requirements for the classical design approach are: 

a. High frequency attenuation In actuator command branch > 15 d6 at 
100 rad/sec and zero angles of attack and sideslip. 

b. Relative stability: gain margins > 6 dB, phase margins >30 deg with 
a goal of 12 dB and 50 degrees. 

c. Acceleration time response: 

(1) 63% time constant of approximately 0.4 sec for a step command of 
acceleration (%) at flight conditions of Interest and at % and 6=0. 

(2) Overshoot! 10% for better sideslip control. 

(3) Steady state error not necessarily zero. 

3. Transfer Functions of Aerodynamic Model and Analysis 

The aerodynamic transfer functions of the uncoupled yaw channel 
autopilot are according to Karadimas, C. thesis [Ref. 41: 

a. Transfer function of yaw angular rate: 

1 

0.15 { s + 1} 

r 0.0748 

= (deg/sec/deg) (IV.B.3- 1 ) 

8y 1 1 

{ s + 1H s + 1} 

4.3 -4.25 

b. Transfer function of Vaw Normal Acceleration: 
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17.76{ 



1 



S+ 1 }{ 



1 



S + 1} 



T\^ 11.48 -11.48 

= (gVdeg) (IV.B.3-2) 

8y 1 1 

{ s ♦ 1} { s ♦ 1} 

4.3 -4.25 

Rearranging the above transfer functions in terms of the variable 
pairs (r,Sv) and (itvySy) respectively and applying the inverse Laplace 

transformation a set of linear differential equations can be obtained. 
Then, using rules of state-space representation, the following equations 
are obtained [Ref. 4]: 



Xt = X 2 - 29.2321 8y (IV.B.3-3) 

X 2 = -16.7224X| - 0 . 0434 X 2 ♦ 0.26875v (IV.B.3-4) 

= Z 2 - 0.080 ISy (IV.B.3-5) 

Z 2 = -1 6.7224 ^ - O.O 435 Z 2 - 150.77248y (IV.B.3-6) 

Y = 3.3574Z, -4Y + 6.434 ISy - 3.3574%^ (IV.B.3-7) 

8Yc= 6.06X^ + O.6O6X2 + 2.04 13Z, + 3.648Y - 1 3.86 1 2 

-2.0413% (IV.B.3-8) 

8Y = 188.48yc - 188.4Sy (IV.B.3-9) 



Using state-space representation, the Equations (IV.B.3-3) through 
(IV.B.3-9) can be modeled in a seventh-order system of the form: 

X = FX ♦ GU 
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with the continuous plant system F and control matrix G as follows: 



Xl' 




0 


1 


0 


0 


0 


0 


-29.2321' 




Xi' 


X2 




-16.722 


-0.0434 


0 


0 


0 


0 


0.26870 




X 2 


Zi 




0 


0 


0 


1 


0 


0 


-0.08010 




Zl 


Z 2 


= 


0 


0 


16.722 


-0.0435 


0 


0 


-150.772 




Z 2 


Y 




0 


0 


3.357 


0 


-4 


0 


6.434 




Y 


Svc 




6.06 


0.606 


2.041 


0 


3.648 


0 


-13.86 




6yc 


l«vj 




0 


0 


0 


0 


0 


188.4 


-188.4 

4 




Sy 

V • / 



■ + 



0 

0 

0 

0 

3.3574 

2.0413 



fnvcl 

I J 



where the state variables are: 

- Xi, X 2 : yaw angular rates 

- Zt » Z 2 : yaw normal acceleration 

- Y : output of acceleration compensator 

- Syc : Input command In the actuator 

- Sy : yaw tall Incidence 



Using a step Input function which represents the "I gee command' 
at ae = 0, the time response plots (Figure 4.4), the Yaw Normal 
Acceleraiion (T\y) (Figure 4.5), and the Yaw Angular Rate (R) and the Yaw 
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Tail Incidence (Sy) (in Figure 4.6) were obtained. The transient response 

satisfies the performance criteria and definitions given at Section (III.B.4}, 
and are identical with the results presented in [Ref. 4). 

In Figure 4.4 the Vaw Normal Acceleration is plotted which has: 

(1) A 0.31 sec time constant. 

(2) A 0.3 sec delay time and 6A% overshoot. 

(3) A 0.018 steady-state error. 

(4) Note the faulty behavior at the start of the response (0 to 0.2 sec). 
Because the system Is a nominimum phase system (i.e., a zero lies 
in the right-half s plane) the transient response starts out in the 
oppsite direction to the input i.e., this is a inherent delay in the 
transient response of the system. 

In Figure 4.5 the Yaw Angular Rate is presented which has: 

(1) A maximum body angular rate Of 0.6 rad/sec achieved at 0.3 second. 

(2) A minor undershoot and 2.2 sec settling time after which the 
angular rate remains constant at 0.01 rad/sec. This is due to the 
fact that the transfer function r/6y (Equation IV.B.3-1) is a type 
zero system and has a finite error to a step input. 

Finally Figure 4.6 shows that a negative angle of deflection of 
-0.14 rad is required in order to achieve the commanded step input 
acceleration. 

The above analysis has revealed that the performance of the system 
confirms to the previous mentioned requirements. 

From the pole-zero plot in Figure 2.7 [Ref.4] we can see that the 
continuous open loop system is stable, since the s-plane poles are: 

Si = -175.010 -jO (IV.B.3-10) 
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Figure 4.4 Yaw Normal Acceleration vs. Time; Uncoupled Yaw Channel 
Autopilot; Continuous Open Loop System; Elliptical Airframe; 
Step Input HY = 1 gee; ag = 0. 
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$2 = -5.7665 +J 10.4804 
S 3 = -5.7665 "J 10.4804 
84 = -0.0217473 +J4.08919 
S 5 = -0.0217473 -J 4.089 19 
Se = -2.9484 +J 2.99489 
S 7 = -2.9484 -J2.99489 



(IV.B.3-1 1 ) 
(IV.B.3-12) 
(IV.B.3-1 3) 
(IV.B.3-1 4) 
(IV.B.3-1 5) 
(IV.B.3-1 6) 



C. UNCOUPLED LINEAR OPEN LOOP ROLL CHANNEL (ORIGINAL) 

1. The Aerodpamic Model, the Airframe Configuration and the 

Control Law 

Generally, the purpose of a roll channel autopilot is to command 
a roll ( 9 ) signal, in order to roll the missile through an angle ( 9 ) 
measured from the vertical axis. 

While the pitch channel acceleration is commanded to produce 
the total magnitude of the guidance acceleration command, the roll 
channel of the B.T.T. autopilot Is commanded to roll the missile so ^s to 
put the prefered maneuver direction in the direction of the guidance 
acceleration command, as is shown in Figure 3.2. 

A block diagram of the uncoupled roll channel is shown in Figure 
4.7, and the roll control law in Figure 4.8 ( both Figures taken from 
[Ref. ID. In these diagrams the same aerodynamic model is used. 

As is shown in Figures 4.7 and 4.8 the control law is commanded 
by a roll angle (9c)i is governed by roll angular rate (P) and roll angle (9) 
and determines the required command ( 61 ^ 0 ) to an actuator. This is 
approximated as a first order lag at 188.4 rad/sec. 
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The circular airframe configuration Is chosen for this channel 
and Is shown in Figure 3.1. The flight conditions are the same as for the . 
Yaw and Pitch channels, i.e. 60000 ft and Mack number 3.95. 

2. Requirements for the Classical Design 
Requirements for the classical design approach are: 

a. High frequency attenuation In Actuator command branch > 15 dB at 
100 rad/sec and zero angle-of-at tack. This requirement limits the 
speed of roll angle response. 

b. Relative stability: gain margins > 30^ with a goal of 12 dB and 50^^, 
I.e. the same as for the yaw channel. 

c. Acceleration time response : 

(1) 63X time constant of Q.5 seconds for a step command of roll 
angle at the flight condition of interest and zero angle of attack. 

(2) Overshoot < lOX. This requirement is conservative and may be 
relaxed depending upon guidance results. - 

(3) Zero steady state roll angle error. If not zero; It will affect 
guidance. 

3. Transfer Funclons of the Aerodynamic Model and Analysis 

The aerodynamic transfer function of the uncoupled roll channel 
autopilot are (Ref. 4). 

a. Transfer functions of the Roll Angular rate: 

P = (^Sd/Ixx)’< 57.3 ^ C| (rad/sec^) 

where: 

^ = 1650 (lbs/ft2) 

5=3. 14159 ft^ 
d = 2 ft 
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(IV.C.3-1) 



Ija =40 slug-ft^ 

C| ^ =0.03 1 
P = 8.0346 6r 

b. Transfer function of Roll Angle: 

= P (rad/sec) (IV.C.3-2) 

Rearranging the above equations along with those given in 
[Ref. 1] for the actuator and using rules of state-space representation, the 
following equations are obtained [Ref. 4] and verified: 



P = 8.034626r (IV.C.3-3) 

=P (IV.C.3-4) 

X = -17.6'f» -8X +17.69o (IV.C.3-5) 

Y = -0.75325P -0.88359 +0.35 165X -5Y -0.403345r 

♦0.883529c (IV.C.3-6) 

X, = -6X, +0.6282291 Sr (IV.C.3-7) 

^0= -0.1 027 IP -0.120489 +0.04795X +14.3182Y 

-14.18184X^ -15Sro (IV.C.3-8) 

Sr = 10795.32Src -188.4Sr (IV.C.3-9) 



Utilizing state-state representation, the above Equations can be 
modeled as a seventh-order system in the format X = FX +GU with the 
continuous plant system F and control matrix G as shown in Table 4.1 and 
the state vector with the corresponding state variables as follows. 
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PLANT SYSTEM AND INPUT MATRICES; UNCOUPLED ORIGINAL 
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p 

9 

X 

Y 

Xi 

Src 



Roll angular rate 
Roll angle 

Output of roll angle compansator 
Output of roll rate compansator 
Ouput of pseudodifferenslator 
Input command in the actuator 
Roll tail incidence 



Executing the interactive implementation of the optimal systems 
control design program (OPTSYSX) using an input step function which 
represents ‘1 radian command* at zero trim angle of attack, the graphical 
time response plots were obtained for the state variables of interest. 

Figure 4.9 shows the Roll Angular Rate (P) (rad/sec vs. time), 
Figure 4.10 the Roll Angle (9) (rad vs. time), and Figure 4.1 1 the Roll Tail 
Incidence (deg. vs. time). A close observation in Figure 4.10 indicates 

that for the commanded input vector the time response plot of the roll 
angle has a 0.55 sec time constant, 5.3% overshoot and a steady-state 
error equal to zero. These results are in accordance with the requirements 
mentioned in the previous section. 

From the pole-zero plot of Figure 2.23 (Ref.4l we can mention 
here that the continuous open loop system is stable, since the s-plane 
poles are: 



s, = -174.785 -jO 
S2 = -9.2510 +j28.4098 



(IV.C.3-10) 

(IV.C.3-11) 



S 7 = -5.19932 +jO 



S 3 = -9.25097 -j28.4098 
84 = -2.46608 ♦J2.71152 
$5 = -2.46608 +j2.7n52 
Se = -8.98209 +J0 



(IV.C.3-12) 

(IV.C.3-13) 

(IV.C.3-14) 

(IV.C.3-15) 

(IV.C.3-16) 



4. Design Aporooch and Analysis of the Modified Model (Unimpaired ) 

To study reconflgurable flight controls the missile model should 
be capable of using the set of four available control surfaces to produce the 
required forces and moments. 

To achieve this each control surface (8), 82, 83, 84) was considered 
as an Independent component. 

The equations were then coupled so that a component of control 
surfaces used alone or In conjunction with another control surface, will 
produce the appropriate moments. 

Taking under consideration the conventional control surface theory 
of Chapter II. for the Bank-to-Tum missiles, and for the circular 
configuration shown In Figure 3.4, one can easily show, for the continuous 
open loop roll autopilot that the Input command In the actuator network 
( 8 pc) can be split up into Independent components, using the following 
equation which computes the deflection Inputs to Equation (IV.C.3-8). Also 
refer to Figures 2.2 and 3.4 which show the control surface position and 
corresponding positive deflections. 






(IV.C.4-1) 
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Replacing the input vector in Equation (IV.C.3-0) and (IV.C.3-9) 
with the above Equation and applying the superposition principle, gives the 
following modified airframe equations: 

= -188.46r + 10795.32 - S 2 - * 6^ / A (IV.C.4-2) 

(6, -62 -63 +64) / 4 = -0.1 027 IP -0.12048f +0.04795X +14.318Y 

-HI8I8X1 -15/4 (61-82-63+64) - 0.140676r 
- 0.12048'Pc (IV.C.4-3) 

In Figure 4.6 the modification of the commanded tail incidence is 
presented. 





^ ACTUATOR 



4 6 



R‘ 



Figure 4.6 (a) Original System; (b) Modified System. 



Rearranging the above equations and using the transfer functions 
of the system for the original uncoupled roll channel, Equations (IV.B.3-3) 



ill 



through (lV.B.3-9), the following Equations were derived and can be 
modeled in a ten order system for the roll channel: 



P, = 8.034628r 


(IV.C4-4) 


=p 


(IV.C.4-5) 


X = -17.69 - 8 X+ 17.69c 


(IV.C.4-6) 


Y = -0.75325P -0.883529 +0.35 165X -5Y 
-0.403345r-+0.883529c 


(IV.C.4-7) 


X| = - 6 X| +0.6829 16r 


(IV.C.4-8) 


= -0.10271P -0.120489 +0.04795X ♦ 14.3182Y 
-14.18184Xt -158^ -0.140628 r +0.120489c 


(IV.C.4-9) 


82 = +0.1 027 IP +0.120489 -0.04795X - 1 4.31 82 Y 
+ 14.18184Xt +158i +0.140628 r -0.120489c 


(IV.C.4-iO) 


83 = +0.1 027 IP +0.120489 -0.04795X - 14.3182Y 
+ 14.18184Xi +158i +0.140628 r -0.120489c 


(IV.C.4-11) 


84 = -0.1 027 IP -0.120489 +0.04795X + 14.3182Y 
-14.18184X, -158i -0.14062 Sr + 0.1 20489c 


(IV.C.4-12) 


8 r = 2698.838, -2698.8382 -2698.8383 +2698.8384 
-188.48r 


(IV.C.4-13) 



We thus have a multiple-input-multiple-output (hlMO) system 
with the following order 

- X : 10x1, the time derivative of X matrix. 

- F : 10x10, the open loop dynamic matrix or plant. 

- X : 10x1, the variables assignments matrix. 

- G : 10x5, the control distribution matrix. 

- U : 5x1, the control input matrix. 
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The above control and output matrices of the Uncoupled Modified 
Unimpaired Roll Channel Autopilot are presented at the Appendix B. 

The Eigenvalues found executing the 0PT5YS Control Program, 
that Is the OET(sI-F) values, show that the open loop system Is stable, 
since the s-plane poles are: 



5, = -185.057 


+ 


JO 


52 = -9.27462 




J29.2679 


53 = -9.27462 


- 


J29.2679 


84 = -2.35821 




J2.70296 


5s = -2.35821 


- 


J2.70296 


5e = -8.89209 


+ 


jo 


57 = -5.18517 




JO 


Se = -15.000 


+ 


JO 


Sg = -15.000 


+ 


JO 


Sifl- 1 5.000 


+ 


JO 



5. Analysis and Comparison 



(roll angular rate) 
(roll angle) 



(Input command In the actuator) 
(roll tall Incidence) 



.. t ne Ori ginal System 



Applying the same step input function as in Section IV.C.3 which 
represents the "1 radian command", and for zero rad deflection command 
for the four control surfaces (i.e, 8 ic 820 ho ^ 4 c = 0 ) the time response 



plots of the following state-variables of interest were obtained and 
presented in the Figure 4.12 through 4.15. 

a. Roll Angle (f or F), Figure 4.12, rad vs. time which has: 

(DA time constant 0.6 seconds 
(2) A 5.3% overshoot. 
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(3) Zero steady-state error and so the guidance will not he effected. 

(4) Delay time 0.48 sec and settling time 2.2 sec when the steady 
state roll angle Is achieved. 

b. Roll Angular Rate (P or pX Figure 4.13, rad/sec vs. time, 
which has: 

(DA maximum body angular rate 1.56 rad/sec at 0.4 second. 

(2) A 4.2% undershoot and a 2.4 sec settling time. 

c. Roll Tail Incidence (6 r or DR), Figure 4.14, degrees vs. time, 



which has: 

(1) Initially and very fast, 0.1 sec, the maximum required deflection 
occurs, that Is 1.17 degrees. Then changes direction to the opposite 
deflection angle In 0.6 sec to -0.32 deg. before stabilizing. 

(2) A 2.6 sec settling time . 

d. Deflections (Sf, 62* h* ^4 of Figure 4.15, rad vs. time, 
which has: 

(1) The actuators Si and S4 have the same sign (positive), as was 
expected according to the sign convention of Figure 2.2, and are 
first deflected to a 0.02 rad angle and then to -0.05 radians after 
0.6 seconds. Finally the steady state zero angle of deflection Is 
achieved at 2.6 seconds. 

(2) The S2 and S3, have the same response sa above but with a negative 
sign. 

As was shown the relative stability for zero-angle-of-attack Is 
satisfactory for the modified unimpaired system. 

A close observation of Figures 4.12 throught 4.15 and the above 
analysis, indicates that the corresponding obtained responses are 
Identical, with those of the original roll channel. Figures 4.9, 4.10 and 
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4.11. 5o, the perTonmance of the designed system meets the pervlously 
mentioned requirements. 

In table 4.2 a summary of the stability margins is presented, 
obtained from the frequency response plots calculated up to 100 rad/sec. 



TABLE 4.2 

PHASE CROSSOVER FREQUENCIES AND GAIN MARGINS; 
UNCOUPLED MODIFIED ROLL CHANNEL(UNIMPAIRED); 
CIRCULAR AIRFRAME 



PHASE CROSSOVER GAIN 
FREQUENCY(RAD/SEC ) MARGIN(DB). 

ROLL ANGLE (f) 10.0461 60.2165 

ROLL ANGULAR RATE (P) 10.0461 60.2165 

ROLL TAIL INCIDENCE (6p) 1 0.046 1 42. 1 1 72 



We procceed now to the analysis of the Impaired control surfaces 
for the modified system. 



D. CONTROL EFFECTIVENESS DUE TO THE CONTROL SURFACES 
IMPAIRMENT AND COMPARISON 

The time response plots for the system are obtained In this section 
for two types of damages. The first type Is the locked or hard over 
actuator at a small angle and the second type Includes the centered or 
Inoperative one, two and three control surfaces. 

Note that the first phase of the design approash for the CBTT 
autopilot was to design each channel Independently with all coupling 
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between channels removed. Sufficient high frequency attennuatlon was 
added for actuator and missile design so that the resulting missile body 
angular rates and control surfaces motion would represent a practical 
missile design. A relationship was established among the relative speeds 
of response of the uncoupled channels In order to meet CBTT requirements. 
The cross-coupling due to control surfaces Impairment between the 
longitudinal and lateral motion of the missile Is Ignored when each control 
channel Is treated Independently. Cross-coupling of the controls Is thus 
considered In Chapter V. 

1. Hard Over One Control Surface 

One control surface (5]), was selected to be deflected or hard 
over, I.e. It remained at a fixed angle while the other three control 
surfaces had a variable deflection angle. 

Using a step Input 1 radian for the commanded roll angle (9c) 

In the original and unimpaired models); an angle of 0.2 rad for the 
damaged surface (6fc); and a zero angle of deflection, for the remaining 
other three (52c; ond54c). The time response plots of the state 
variables of Interest were obtained and analyzed with the performance 
criteria as defined In section (III.B.4). 

a. Roll Angle (9 or F), Figure 4.16, (radians vs. time) has: 

(1) A time constant 0.6 sec and 0.3 sec delay time . 

(2) A constant 0.4 rad steady state error. 

(3) A 5.8* overshoot at 1.2 sec. 

(4) The roll angle Is stabilized at an Incorrect angle of 1.4 rad instead 
of the commanded 1 radian due to the Impairment of 6^ surface. 
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b. Roll Angular Rate ( P ), Figure 4.17 (radians/sec vs. time) 



which has: 

(1) Initially and due to the fact that S) Is locked at a positive angle of 
0.2 rad, to achieve the the roll command, the angular rate Increases 
very fast the first tenth of the second to a 2.1 rad/sec rate. Then, a 
pronounced oscillation occurs until the 0.4 sec when the rate begins 
to decrease crossing the zero rate at 1.2 sec corresponding to the 
max roll angle (1.5 rad). 

(2) The Initial Increase of the rate Is due to the Induced rolling moment 
caused by the deflection of the damaged control surface while the 
other three are at their steady position. 

(3) A 6.6% undershoot at 1.5 sec and settling time 2.5 sec. 

c. Roll Tail Incidence (6q or DR), Figure 4.18 (deg. vs. time) 



which has: 

(1) At time equal to 0.1 sec the actuator Is commanded to an angle 4.75 
degrees. This Is is due to the fact that 6^ Is hard over at the 0.2 
rad angle while the other three remain have a zero Initial condition. 
Then the tall Incidence Is stabilized after a 2 second period during 
which a fluctuation occurs. 

(2) A zero steady-state error and 2.2 sec settling time. 

d. Control Surface Deflections (Sf, Sz, $3,64 or D's), Figure 
4.19, rad vs. time, which has: 

(1) The actuator Sf, Is deflected at the angle of 0.2 rad, from zero sec 
up to 3 seconds. 

(2) The commanded Input vector now causes the two surfaces Sz and S3 
to turn In a positive angle of 0.13 rad. Instead of a negative angle 
as when no Impairment existed. 

(3) A 1.6 seconds settling time to achieve this angle. 
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(4) The fourth surface 64 Is commanded to turn at an angle of equal 
amount but with an opposite sign, for the same reasons. 

A comparison between Figure 4.19 and 4.15 of the unimpaired 
system, shows that to achieve the roll response the impaired sustem 
requires more control surfaces deflection and more time. 

2. Damage One Control Surface 

The 61 control surface was selected to simulate the Inoperative 
actuator. The corresponding state variable 5i was therefore set equal to 
zero. The step input of one radian remained the same. The initial condition 
on the control surface was set at zero. This type of damage could also 
represent the dedtruction of the control surface. 

For the above case the time response plots of the state variables 
of Interest were obtained and analyzed. Later in Section 5, a comparison 
between the different type of damages and the unimpaired system is 
performed. 

a. Roll Angle ( f or F ), Figure 4.20, rad vs time, which has: 

(1) A 0.6 second time constant. 

(2) An overshoot equal to 5.6% and 

(3) Zero Steady state error. 

(4) Delay time 0.5 sec and 2.1 sec settling time at which the the roll 
response is achieved. 

b. Roll Angular Rate ( P ), Figure 4.21, radian/sec vs. time, 
which has: 

(1) A maximum body angular rate 1.58 rad/sec at 0.4 sec. 

(2) A 6.3% undershoot. 
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(3) A 2.5 sec settling time at which the zero steady state value Is 
achived. 

c. Roll Tail Incidence ( 6 ^ or DR ), Figure 4.22, deg. vs. time, 



which has: 

(1) At the very beginning of the response, 0.1 sec, the maximum 
deflection occurs, that is 1.1 degrees, to change then to the 
opposite deflection angle in 0.7 sec to -0.32 deg. 

(2) A zero steady state error and settling time 2.6 sec. 

d. Control surfaces deflection or DR), Figure 4.23, deg. vs. 

time which has: 

(1) The surface as damaged remains constantly at zero angle of 
deflection. 

(2) The actuators 83 and S 3 are deflected at a negative angle 0.026 rad 
and are stabilized after 2.4 seconds. 

(3) The actuator 84 Is deflected at an equal but opposite angle as It was 
expected according to the selected sign convention. 

A comparison between the deflection angle in Figures 4.15, 4.19 
and 4.23 Indicates that for the same step Input: 

(1) When one control surface Is Inoperative or shot away the system 
commands the undamaged surfaces at larger angles (In comparison 
with the unimpaired system). 

(2) When one control surface Is inoperative or shot away the system 
commands the undamaged surfaces at smaller angles than obtained 
with the hard over type of damage. 

(3) The required time that the undamaged control surfaces need to 
achieve the zero steady state value Is the same (In comparison with 
the unimpaired system). 
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3. Damage Two Control Surfaces 

In this case two control surfaces, on the same side of the 
missile were considered damaged (6) and 62). The two state variables Sy 
and 82 were therefore set equal to zero. The step Input of one radian 
remained the same. The Initial condition of the control surfaces was set 
at zero. Then the time response plots of the state variables were obtained 
and analyzed. 

a. Roll Angle (9 or F), Figure 4.24, rad vs. time, which has: 

(1) Time constant 0.59 seconds 

(2) An overshoot 6.8% at 1.25 sec and an undershoot of 0.003% at 2.2 
seconds. 

(3) A zero steady-state error and a 0.5 sec delay time. 

(4) The commanded roll angle Is achieved after 2.8 seconds. 

b. Roll Angular Rate (P), Figure 4.25, radian/sec vs. time, 
which has: 

(1) The peak body angular rate 1.61 rad/sec Is achieved at 0.41 sec. 

(2) An 8.2% undershoot at 1.6 sec and 

(3) A 2.6 sec settling time. 

c. Roll Tail Incidence ( 8c or DR ), Figure 4.26, deg. vs. time, 
which has: 

(1) The peak angle Is 1.01 degrees at 0.1 sec and then In 0.6 seconds a 
negative angle (-0.35 deg) is commanded. 

(2) A second 6.3% overshoot at 1.8 sec and 

(3) Zero steady state error (settling time 2.6 sec). 
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d. Control surface deflection ( 6 i, 62 * ^ 3 * ^4 or D's), Figure 4.27, 
rad vs. time, which has: 

( 1 ) The two Inoperative or shot away control surfaces 6 | and $2 remain 
always at zero angle of deflection. 

(2) The two remaining actuators S 3 and S 4 moving at opposite 
directions, as expected, with a peak value of ±0.35 radians. 

(3) The Impairment of Sf and 62 has produced a more overshoot lags 
system as evidenced In Figure 4.27. 

(4) The convergent and oscillatory behavior of the two unimpaired 
control surfaces which attempt to "fill in" for the damaged surfaces 
is shown In the same Figure. Also note that for this system with 
two Impaired control surfaces more angle of deflection Is required 
(Incomparison with the system having one actuator shot away). 

4. Damage Three Control Surfaces 

In the last case three of the control surfaces 62 2 nd §3 were 
damaged or Inoperative. 

With a step input of one radian for the commanded roll angle ('Pc) 
as in the unimpaired and other type damaged models, the plots of the state 
variables of interest were obtained and analyzed. 

Figures 4.28 through 4.31 Illustrate the response of the output 
variables. Figure 3.30 shows that the commanded one radian roll angle can 
still be achieved In the absence of three control surfaces (S^, S 2 , 63 ). 
However, some oscllatlon Is evident In the response. The loss of the 
surfaces has produced a more lightly damped system as evidenced in 
Figures 4.28 through 4.31. The oscillatory and convergent behavior of the 
roll rate and the undamaged surface (S 4 ), shown in Figures 4.29 and 4.31, 
Illustrate how the surface attempt to "fill In" for the absence of the 
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Inoperative control surfaces. Specifically the response of the roll angle 
(Figure 4.28) has: 

(1) A 0.55 sec time constant. 

(2) A max. 10.5% overshoot at 1.2 sec and a 0.5% undersoot at 2.2 sec. 

(3) The commanded roll angle is achieved after 2.8 sec. 

(4) A 0.51 sec delay time . 

A comparison with the response of the corresponding deflection of 
the unimpaired system (Figure 4.15) and the other inoperative cases. Figures 
4.19, 4.23 and 4.31, indicates that as the number of the inoperative control 
surfaces increases, the deflection of the undamaged actuator(s) also 
increases and more time is required to acheive th steady-state value. 

5. Comparison with the Unimpaired 5ustem 

In this section a performance comparison of each type of 
impairment with the unimpaired system, was performed and an analysis of 
their corresponding time plots is presented. 

For the above purposes, a computer program was written 
(Appendix 0) to create a data file, for the state variables of interest, to 
be used as input data in the use of the EASYPLOT program of the NP5 
computer. 

For the decoupled roll channel the first plots are for the roll 
angle and are presented in Figures 4.32 and 4.33 

In Figure 4.32 which compares the two types of impairments 
with the original (unimpaired) system it is quite obvious that the worst 
case is the hard over case because when the autopilot is commanded at a 
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roll angle (fc) equal to 1 radian, the missile turns at a roll angle (9) of 
1.4 rad Instead of 1 rad as for the unimpaired autopilot. 

Other characteristic point is that for the hard over case, the roll 
angle reaches its final value of 1.4 rad (80.2^) in approximately 2 seconds 
with a 5.8% overshoot, without any steady state error. 

Figure 4.33 compares the roll angle of the unimpaired system 
with the three cases of control surfaces impairments. The curves include 
the following four cases: 

(1) Original system 

(2) Inoperative one control surface 

(3) Inoperative two control surfaces 

(4) Inoperative three control surfaces 

This Figure shows that as the number of the damaged control 
surfaces increases the overshoot and the settling time icrease. 

Figures 4.34 and 4.35 compare the roll angular rate of the 
unimpaired system with four cases of control surfaces impairments. The 
curves include the following five cases: 

(1) Original system 

(2) Hard over one control surface 

(3) Inoperative one control surface 

(4) Inoperative two control surfaces 

(5) Inoperative three control surfaces 

In Figure 4.34 it is shown that the hard over case induce greater 
excursion in the roll rate with a peak value of 2.2 rad/sec at 0.38 sec 
being achieved. Due to the fact that control surface is deflected at the 
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0.2 rad the Initial rolling moment in the first tenth of the second 
increases from zero to 2.1 rad/sec. The other control surfaces have zero 
as an initial condition. 

Figure 4.35 compares the unimpaired system with three cases of 
control surfaces impairments (inoperative one, two and three actuators). 
In this is shown that the inoperative cases are all quite similar except for 
the last case where only one control surface is available. In this last case 
a larger delay time, overshoot and undershoot are evident. 

Finally in Figures 4.36 and 4.37 the roll tail incidence is 
compared and include the following cases: 

(1) Original system 

(2) Hard over one control surface 

(3) Inoperative one control surface 

(4) Inoperative two control surfaces 

(5) Inoperative three control surfaces 

Figure 4.36 shows that the required roll incidence for the hard 
over case is larger in comparison with the unimpaired system and with 
the other three cases of control surfaces impairment. In Figure 4.37 the 
maximum peak overshoot decreases with the number of control surfaces 
that are impaired. 



133 



F.ROLL ANGLE-RAD 




134 




P.ROLL ANGULAR RATE-RAD/SEC 

- 0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 




135 



DR, ROLL TAIL INCIDENCE- DEG. 

0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 



o 




Figure 4.18 Roll Tail Incidence vs. Time; Uncoupled Modified Roll Channel 
with Hard Over ; Continuous Open Loop System; 

Step Inputs: 5i=0.2 rad; 52 , 63 , 64 = 0 ; fc=1 rad and a^ = 0. 
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Figure 4.19 Control Surfaces deflections vs.time; Uncoupled Modified Roll 
Channel with Hard Over ; Continuous Open Loop System; 
Step Inputs: 5t=0.2 rad; 62 , 63 , 64 = 0 ; ^nd ag = 0. 
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Figure 4.21 Roll Angular Rate vs. Time; Uncoupled Modified Roll Channel 
with Centered; Continuous Open Loop System; 

Step Input: 9c=l r^d and Sic»fi 2 c»^ 3 c*^ 4 c"^l ~ 
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Figure 4.36 Comparison of the Roll Tail Incidence between the 

Unimpaired System and the Four Cases of Impairment. 
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V. COUPLED LINEAR PITCH AND ROLL CHANNELS 



A. INTRODUCTION 

Separate longitudinal and lateral autopilots were designed in Chapters 
III and IV, for the Modified Linear Open Loop system. 

In this Chapter the longitudinal and roll autopilots were combined, 
resulting in a new autopilot capable of commanding pitch, roll or combined 
maneuvers. 

The time responses to an input vector are identical to that obtained 
when the same input is put into the separate pitch and roll autopilots of the 
unimpaired system. 

The response plots of this system are then used as a basis to 
investigate the effect of the damaged control surfaces on the combined 
system. In this section the hard over type of damage for the combined 
system is investigated. Corrective actions are then developed and presented 
and the robustness of the complete system is finally analyzed. 

B. DESIGN APPROACH, ANALYSIS AND PERFORMANSE OF THE SYSTEM 

For the design of the combined system of the circular airframe 
configuration, the requirements and assumptions mentioned in the design of 
the separate autopilots are used. To achieve the desired maneuver plane 
acceleration response for the unimpaired system, both pitch and roll 
uncoupled autopilots were designed to have the desired maneuver plane 
speed of response. In paticular, the goal for maneuver plane acceleration 
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response was a 0.5 second time constant for the flight condition of Interest 
(I.e., 60 Kft altitude, Mach 3.95. in the next phase of the linear design 
technique the longitudinal and lateral autopilot channels are coupled using 
the aerodynamic model of Chapter III. Both channels are cross-coupled 
through the control surfaces (I.e., a command to the control surface of one 
channel will affect the parameters of the other and vice versa). The pitch 
and roll state Equations (III.B.5-5) through (III.B.5-17) and (IV.C.4-4) 
through (IV.C.4-13) are comdined to yield a twenty-three order system 
which Is presented In Appendix C. 

The Open Loop dynamic matrix Is a 23 x 23 matrix. This composite 
system Is formed by two block diagonal matrices as Is stated In Equation 
(V.B-1): 



F| (13x13) I 0 



F = 



I 

0 I Fz (10x10) 



(V.B.1-1) 



That is, the upper block composed of the 13x13 F matrix from the 
pitch channel and the lower block composed of the 10x10 F matrix from the 
roll channel. 

The control distribution matrix G, folowing the same pattern, is given 
by a 23x10 matrix, of the form and order: 



G| (13x5) 


0 


0 


62 (10x5) 



(V.B.1-2) 
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Finally the state vector X and the control input vector U matrices are 
of the form and order. 

X= I23x 1 1 (V.B.1-3) 

U= I lOx 1 1 (V.B.1-4) 

Because of this form, any input to the stated combined system, has a 
time response identical to that resulting from the same input into the pitch 
and roll channels separately. An example of this is shown in Figure 5.1 
where the same input vector used for the vertical translation maneuver 
longitudinally and the rolling translation maneuver laterally are combined 
and serve as inputs to the combined system resulting in a response 
identical to those achieved separately in Figures 3.9 through 3.12 and 4.12 
through 4.15. 

The combined autopilot designed in this thesis and the simulations 
made to test the different step input commands are accomplished using the 
interactive, multivariable design program called OPTSYS. 

in addition, in order to compare different control surface impairment 
configurations a computer program was developed (Appendix 0). This 
program creates data files to be used as input files at the EASYPLOT 
program by changing the format of the OPTPLOT data into a fixed format 
data file. This Format is 80 characters in length with space between the 
values. 

Three different command inputs are examined for the Modified 
Unimpaired system and are discussed below. 
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Test 1 : For this lest, both a pitch normal acceleration and a roll angle 
were commended (i.e., = 1 gee, ^ ^ic» ^ 2 o ^ 3 c» ^ 4 c = 

In Figure 5.1 the time response plots of the normal acceleration and the roll 
angle are presented. As was expected, both curves are identical to the 
corresponding plots obtained when separate autopilots were used. 

Figures 5.2 through 5.4 illustrate the performance of the other state 
variables. Figure 5.4 shows the control surface deflections for only two 
control surfaces from each channel. The remaining four coincide or are the 
mirror image of the obtained results. In the same Figure the configuration of 
the control surfaces for each channel is provided and the distinction is also 
provided in the legend. 

Test 2: For this test, only a pitch command in acceleration was applied 
( i.e., = 1 gee), while , 9c = ^ rad and 620 ^30 ^ 4 c = ^ rad. Figures 

5.5 through 5.8 show the performance of the system. Observe that no roll 
angle or any roll rate is present, as was expected, while a downwards 
maneuver is performed according to the command input vector ( downward 
maneuver is the positive one). Figure 5.8 shows the response of the control 
surfaces for both channels. The control surfaces of the longitudinal channel 
are deflected at a ±0.025 radian angle and of the lateral channel remained at 

the initial zero position while only a downward maneuver was commanded. 
Test 3 : In the last test, only a roll command is applied ( i.e., 9c = 1 

), while N^c = 0 gee and 620 ^30 ^ 4 c = ^ rad. Figures 5.9 through 5.12 

Illustrate the effectiveness of input vector on the state variables of 
Interest. In Figure 5.12 only the two control surfaces 5^ and 82 of the roll 
channel are plotted while the corresponding surfaces in the pitch channel 
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remain at zero. The other four coincide with the obtained results. The 
control surfaces of the roll channel are deflected at a ±0.022 radian angle 
and of the longitudinal channel remained at the Initial zero position while 
only a rolling maneuver was commanded. 

it should be mentioned here that according to the pole-zero plot for the 
designed combined system, the combined continuous open loop system Is 
also stable, since the s-plane poles are: 



51 = -159.724 ♦ j 18.9922 (V.B-5) 

52 = - 1 59.724 - j 1 8.9922 (V.B-6) 

Sj = -8.29043 + jO.0591 1 (V.B-7) 

S 4 = -8.29043 - jO.0591 1 (pitch angular rate) (V.B-8) 

Ss = -3.75933 + J2.51441 (V.B-9) 

Se = -3.75933 - J2.51441 (V.B-10) 

S7 = -0.14393 +jO (V.B-11) 

Sg = -0.00005 + jO (pitch tall Incidence) (V.B-12) 

Sg = -0.07460 + jO.66599 (V.B-13) 

Sio= -0.07460 - jO.66599 (pitch normal acceleration) (V.B-14) 

S|i= -0.14300 +jO (V.B-15) 

512 = -0.14300 +jO (V.B-16) 

513 = -0.14300 + jO (V.B-17) 

51 4 = -0.01850 + jO (roll angular rate) (V.B-18) 

Sis = -9.27462 + J29.2679 (roll angle) (V.B-19) 

Sie = -9.27462 - J29.2679 (V.B-20) 

Si7 = -2.35821 +J2.70296 (V.B-21) 
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S2J ~ ~ 1 5.0000 jO 
$22 = -15.0000 + jO 
$23 - "15.0000 jO 



$t8 = -2.35821 -J2.70296 
$,9 = -8.89209 + jO 

$20 = -5.18517 + jO (roll tail incidence) 



(V.B-22) 

(V.B-23) 

(V.B-24) 

(V.B-25) 

(V.B-26) 

(V.B-27) 



C. THE EFFECTS DUE TO ACTUATOR IMPAIRMENT 

In this section the effects on the system due to a damaged actuator 
was investigated. One control surface (Sf) was considered damaged. It was 
locked at a fixed angle while the other three control surfaces had a variable 
deflection angle. 

Two different control inputs were examined for this type of 
impairment and are discussed below. For each test, corrective actions were 
then developed and presented. 

Test 1 . Command Input: 

(1) for both channels locked at 0.2 rad 

(2) 52, 63, 54 for both channels have variable deflections 

(3) Nzc = 1 gee 

(4) 'pc = 0 rad 

Using the above control input the time response plots of the system 
were obtained, analyzed and presented in Figures 5.13 through 5.16. in 
Figure 5.13 the pitch normal acceleration and the roll angle are presented. 
Commanding a downward maneuver (l.e., Nz=l gee and 'p = 0) the missile 
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body the first tenth of the second moves downward and then upward with its 
pitch axis directed at the velocity vector until it reaches the desired 
maneuver level. This is due to the Induced pitching moment effect of the 
locked control surface and the undamaged control surfaces (62> S?. S4) 
which try to counteract the impairment of AlUhough no roll motion is 
required to maintain cordination for this maneuver, as is shown in Figure 
5.13 a roll angle (0.4 rad) is present. In Figure 5.14 the pitch and roll 
angular rates are also presented in which the 1.9 rad/sec roll rate Is 
achieved. In Figure 5.16 two control control surfaces from each channel are 
presented. The other two are the mirror images of the obtained results. The 

6 1 (of pitch and roll channels) respond properly to the command input. The 

62 for the pitch channel is deflected 0.174 radians and for the roll channel 
0.130 radians. Both are stabilized after 1.4 seconds. 

Corrective actions : The roll effect of the damaged control surface on 
the performance of the combined autopilot was examined here and 
corrective actions to overcome this problem were studied and presented. 

Figure 5.13 shows the roll angle (9) reaching Its final value of 0.40 rad 
(230) in approximately 1.4 seconds with a 7 % overshoot. The first try to 
overcome the roll effect was to command the opposite surface (83) in an 
equal angle of deflection 0.2 radian. The obtained time response of the pitch 
normal acceleration and the roll angle are presented in Figure 5.17. For 
comparison. Figure 5.17 presents three configurations. It shows that the 
produced roll angle due to the 81 control surface impairment is eliminated. 
But, there Is a penalty In the pitch normal acceleration which reaches 
approximately 4 gee and a oscillatory behavior Is evident before the steady 
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state value of 1 gee command being achieved. So, another test was 
performed by commanding the 63 surface at an opposite angle of deflection 
(-0.2 rad). By this action, shown In Figure 5.18, the normal acceleration 
responds very well but the roll angle Is Increased at 1.18 radians (68°) after 
1.75 seconds, as might be expected. The next test was performed using the 
second control surface (S2). In Figure 5.19 the time response of the same 
variables of the system are obtained and presented. The control surface 62 
Is commanded at the same angle of deflection I.e., 0.2 radians. As Is 
lllusrtated In Figure 5.19 commanding a downward maneuver the missile 
body reaches the desired maneuver level without any roll motion (I.e., Nj =1 

gee and 9 = 0). 

At the end of this Chapter a transition reconfiguration strategy Is 
proposed addressing the topics of failure detection and Isolation, and 
Implementation. That Is, one. pssible strategy for transitioning from the 
normal control laws for a missile with all control surfaces operative, to 
one set of reconfigured control laws for a missile with a failed surface. 

Test 2 . Command Input: 

(1) S) for both channels locked at 0.2 rad. 

(2) S2> ^4 channels have variable deflections. 

(3) Nzc=Ogee. 

(4) 9c = I rad. 

Using the above step Input vector the time response plots of the system 
were obtained and presented In Figures 5.20 through 5.23. 

In Figure 5.20 the time response of the pitch normal acceleration and 
the roll angle of the combined system are Illustrated. Due to the fact that 
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the two channel autopilots were cross-coupled through the four control 
surfaces {8y, S 2 , 8^, 64 ), a damage in any actuatoKs) will affect the 
performance of both channels. The locked control surface 8^ at the one rad 
angle affects the combined autopilot. A constant incorrect level was 
achived in the response of the roll angle (i.e., = 1.4 rad instead of the 

commanded one rad), in addition, due to the Induced pitching moment effect 
by the damaged control surface, a down and upward maneuver was 
commanded by the pitch channel during the first second of the response of 
the normal acceleration. 

A change of the other state variables was obvious as is shown in 
Figures 5.21 through 5.23. Figure 5.23, which represents the time response 
of two control surfaces from each channel, i.e., 6 ^ and 62 , indicates the 
commanded (simulated) failure of 0.2 radian at 8 ^ and the deflection angle 
of 0.13 radians (7.5°) for the other three control surfaces. The control 
surfaces not presented in this plot are the mirror images of the obtained 
ones. 

Corrective actions: The roll and pitch effects caused by the damaged 
control surface 61 on the combined autopilot performance was examined and 
corrective actions to overcome this failure were studied and presented. 

In Figure 5.20 the undesirable rolling and pitching effect caused by 
the damaged St control surface is obvious. The missile rolls about its 
velocity vector with a 1.4 radians roll angle ( 8 O 0 ) (instead of the 
commanded 1 radian). Also observe the effect on the pitch normal 
acceleration which should be constant at zero since no down or upward 
maneuver was commanded. 
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To overcome these failures several runs of the Controls Program 
0PT5YS were performed to find an appropriate corrective action for the 
control system. The most Interesting cases are Illustrated In Figures 5.24 
and 5.25. 

In Figure 5.24 the roll angle (9) and the pitch normal acceleration (n^) 
were obtained and presented for three cases to compare their responses. A 
comparison was performed between the original system, the impaired 
system with the effector locked at an angle of 0.2 radians, and the 
system with the 62 control surface commanded at an angle equal but 
opposite (-0.2) to the defective surface. As Is Illustrated in Figure 5.24, if 
such an Input was commanded, the missile body rolls even more at a 2.2 
radians (126®) angle as might be expected. Also, the pitch normal 
acceleration increases and achieves the large value of -4.2 gee's before 
stabilizing after 1.4 seconds. 

In Figure 5.25 the roll angle (f) and the pitch normal acceleration (nz) 
were obtained and presented but now the second control surface 62 was 
commanded to be deflected at the same angle as the (0.2 radians). As is 
illustrated in Figure 5.25, both effects, on the roll angle and the pitch 
normal acceleration, were eliminated and their time response plots 
Indicates that the missile body performed according to the command input. 

in comparing output time responses, as was done by these two tests. It 
Is apparent that the Impaired missile was able to track the unimpaired 
missile as long as there was another control surface to provide the 
appropriate corrective action and control power of the failed surface. 
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In the case where one or more control surfaces are inoperative or shot 
away and are not considered here, similar effects are expected. This Is due 
to the fact that as far as all the control surfaces are operative, are 
commanded by the system such that one will counteract the induced moment 
effect of the other. But when one control surface Is Inoperative or shot 
away undesirable roll, pitch or combined effects will be present. 

D. PROPOSED RECONFIGURATION STRATEGY 

1. int roduct i on 

This section presents one possible strategy for transitioning from 
the control laws for a missile with all surfaces operative, to a set of 
reconfigured control laws for an missile with a failed surface. The 
discussion begins with a definition of reconfiguration and why It Is 
Important. A transition stategy Is then proposed, addressing the topics of 
failure detection. Isolation, and Implementation. 

2. Reconfiguration 

Reconfiguration Is a technique which shows promise in restoring 
missile stability and a degraded level of controllability in the event a 
control surface becomes Inoperative. This technique Implies that the control 
laws governing missile motion compensate for the failed surface by using 
the remaining surfaces to maintain control. The control laws of Chapter ill, 
IV and the corrective actions of Chapter V provide the Independent control 
of each surface necessary to Implement reconfiguration which deals with 
the coupled aerodynamics. 
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3. Transition Straleq u 

The development of reconflgurable control laws, such as those 
developed In this thesis, must Include a method of failure Identification and 
Implementation of a revised control law. For a given failure It thus 
necessary to determine when the new laws should be used and how to 
switch to them from from normal control laws. A failure detection and 
Isolation scheme must be developed to determine when and which surface 
has failed. The proposed failure detection and isolation system would use a 
set of redundant sensors and an estimator, to determine If the actuators are 
still operative. Should an actuator fall to move a control surface, the 
sensors notify a control surface monitor. The monitor validates the fall 
signal and determines which control surface Is Inoperative. 

The control surface monitor is the heart of the reconflgurable flight 
control system. When notified of a failure, the monitor Informs the flight 
control computers of the problem and provides the necessary Information 
concerning which set of new control laws (corrective actions) should be 
used to calculate the appropriate surface deflections. The off-line flight 
control computer is programmed with the new control laws and then signals 
the on-line system that It Is ready to assume control of the missile. Thus 
the autopilot Is again ready to control a stable missile with some degree of 
control capability. 



179 



N.ACCELER-G’S AND ROLL ANGLE-RAD 

1.50 - 1.25 - 1.00 - 0.75 - 0.50 - 0.25 0.00 0.25 0.50 0.75 1.00 1.25 




180 




ANGULAR RATE-RAD/SEC 

0.4 - 0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 




161 



TAIL INCIDENCE- RAD - DEG. 

- 2.5 - 2.0 - 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 




182 



DEFLECTIONS-RAD 

0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 



A n fmn n nnr! nR^i R PHT i ■ r: -' '7Mr! n fTnnn 




□ 


LEGEND 

Dl-PITCH-RAD 


0 


Dl-ROLL-RAD 


A 


D2-R0LL-RAD 


-f 


D2-PITCH-RAD 





_ ^ 1 1 1 1 ^ 

0.0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 

TIME-SEC 

Figure 5.16 Hard Over 5^ Type of Impairment: 

Pitch and Roll Deflections vs. Time; Combined Pitch 

and Roll Channels; Continuous Open Loop System; Step Inputs: 

5ic=0.2rad, Sjc 53 ^, 64 c = ^ and N^c = 1 gee, 9c = 0 rad; ag = 0 . 
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Figure 5.17 First Corrective Action and Comparison of the Pitch 

Normal Acceleration (Nz) and the Roll Angle ( 9 ), between 
Three Different Effectors Configuration. 
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Figure 5.18 Second Corrective Action and Comparison of the Pitch 

Normal Acceleration (Nz) and the Roll Angle (<f), between 
Three Different Effectors Configuration. 
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Figure 5.19 Final Corrective Action and Comparison of the Pitch 

Normal Acceleration (Nz) and the Roll Angle (9)> between 
Three Different Effectors Configuration. 
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Figure 5.23 Hard Over Type of Impairment; 

Pitch and Roll Deflections vs. Time; Combined Pitch 

and Roll Channels; Continuous Open Loop System; Step Inputs: 

S|c=0.2rad, 62 c, ho S 4 c = ^ = 0 gee, 9c = 1 rad; ag = 0. 
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Figure 5.24 First Corrective Action and Comparison of the Pitch 

Normal Acceleration (Nz) and the Roll Angle (»P), between 
Three Different Effectors Configuration. 
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Figure 5.25 Second Corrective Action and Comparison of the Pitch 

Normal Acceleration (Nz) and the Roll Angle (f), between 
Three Different Effectors Configuration. 
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E. ROBUSTNESS OF THE UNIMPAIRED SYSTEM 
1. Introduction 

Up to this point the continuous open loop system was investigated 
and it was found to De stable, as is indicated by Equations (V.B-5) through 
(V.B-27). In this section the concept of feedback is introduced and the 
robustness of the system is examined. 

Design criteria for feedback control systems can be lumped into 
three broad categories: performance, stability, and robustness. 

Performance criteria are generally specified as a requirement for the 
control system to follow commands or reject disturbances. The most 
commonly used stability criteria are stability margins which generally 
refer to the amount of gain or phase variations from a design condition a 
system employing feedback control can experience before becoming 
unstable. 

in early feedback control applications, reducing sensitivity was 
the only control design criteria, particularly reducing sensitivity to noise. 
Over the last several years, sensitivity Issues have attracted new attention 
and have generated a new requirement, robustness. A control system Is 
robust (Insensitive) If performance Is maintained In the presence of 
variations in the plant from the design model. These variations can be 
modeling approximations, modeling uncertainties, failures of control 
elements (sensors, actuators), noise, non-linearities, etc. 

These criteria are not mutually exclusive. Good stability 
characteristics In well designed systems are usually accompanied by good 
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performance. Robustness must be considered relative to the vacations from 
the design model that the control system will experience. 

2. The Benefits of Feedback and the Singular Value Design Technique 



d 




Figure 5.26 Standard Feedback Configuration 



We deal with the standard feedback configuration illustrated in 
Figure 5.26. It consists of the interconnected plant (G) and controller (K) 
forced by commands (r), measurement noise (n), and disturbances (d). 

The noise in this cose is assumed to be zero. The disturbances ore 
assumed to be reflected to the measured output (y). 

All signals are multivariable, in general, and both nominal 
mathematical models for 6 and K are finite dimensional linear time 
invariant systems with transfer functions matrices G(s) and K(s). Then it is 
well known that the configuration, if it is stable, has the following major 
properties: 

a. Input-output behavior 

y = G K (I + G K)-Ur - n) + (I + G K)-« d (V.E.2-1) 
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(V.E.2-2) 



e = r-y 

= (I + GK)-' (r-d) + GK(I + 6K)-' n 

b. System sensitivity 

AH cl = (I + K)- 1 AH 0, ( V.E.2-3) 

In Equation (V.E.2-3), AH d and AH d denote changes in the closed 

loop system and changes in a nominally equivalent open loop system, 
respectively, caused by changes in the plant G, i.e., G* = G + AG. 

Equations (V.E.2-1) through (V.E.2-3) summarize the fundamental 
benefits and design objectives Inherent In feedback loops. Specifically, 
Equation (V.E.2-2) shows that the loop's errors In the presence of commands 
and disturbances can be made -small" by making the sensitivity operator, 
(I + G K)"*, 'smair. Equation (V.E.2-3) shows that loop sensitivity Is 
Improved under these same conditions, provided AG is not too large 
from G. 

With the increased Interest in MIMO systems numerous methods of 
design have been employed to obtain suitable system performance and 
robustness with varying degrees of success. One primary method of design 
Is to keep the plant as decoupled as possible throughout the design so that 
each Individual element may be controlled Independently and designed 
essentially as a single loop system. In another procedure the multiloop 
system Is modified Into a system that has diagonal elements that are much 
larger than any off-diagonal elements. This diagonally dominant system Is 
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then In a form where conventional Nyquist type techniques can be employed 
In the analysis. A third common tlltIO design method Is that of the Linear 
Quadratic (LQ) method. This method uses a quadratic cost functional and 
optimization principles to allow the designer to design for various 
performance levels by adjusting the matrix weighting terms used In the 
cost function. The major difficulty with all of the above methods Is that 
they are not necessarily robust. This Is especially true for cross-coupling 
terms between loops. 

As was mentioned before, a simple Interpretation of robustness Is 
the ability of the system to tolerate design pertubatlons. These 
perturbations could be In the form of actuator failures, plant parameter 
uncertainty, unmodeled dynamic or nonlinear terms, or any one of many 
other perturbations to the nominal design of the system. 

The analysis of the robustness of the system was obtained using 
the Pole Placement and Robustness Design Program (POPLAR) developed by 
Gordon, V.C. iRef. 51. 

Utilizing the above program and using the data file shown In 
Appendix E for the state-feedback autopilot, the unimpaired system was 
found to be robust. The standard full state-feedback design was carried out 
using the N.P.S. program OPTSYS. 

Figures 5.27 and 5.28 show the ninimum Additive Output Singular 
Value (SVADMO) vs. frequency, and the Minimum Additive Input Singular 
Value (MIN.ADD.IN.SV) vs. frequency respectively for the state-feedback 
designed autopilot. It Is noted that for 0.1 to 100 db frequency the output 
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F*6 Singular value Is too low ( 0.03 rad ) but the Input transfer function 
6*F Input minimum singular value is above the 0.5 rad region. 

Therefore the closed loop unimpaired system is robust (at least in 
the input). Further analysis would be needed if robustness were required in 
both the input and output case simultaneously. 

In the case of state-feedback autopilot and for a singular value of 
0.5 rad the pole placement and robustness design routine places the poles as 
follows: 
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Figure 5.27 SVADMO vs. Frequency; Copied Pitch and Roll Channel 
Autopilot; State-Feedback Design; Continuous Closed 
Loop Unimpaired System; Circular Alrfram. 
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Figure 5.28 MIN ADD IN SV vs. Frequency; Copied Pitch and Roll Channel 
Autopilot; State-Feedback Design; Continuous Closed 
Loop Unimpaired System; Circular Airfram. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The purpose of this thesis was to examine the control deficiecies in a 
Bank-to-Tum (B.T.T.) cruise missile, due to specific impairments at its 
control surfaces during the intermediate phase of its mission. The following 
are the principal conclusions based on this work. 

(1) The original and modified classical designed autopilots which involve 
separate inputs for each control surface were proved to have 
identical performance for the longitudinal and lateral channels. 

(2) The investigation of specific type of damages at the B.T.T. missile s 
control surfaces, showed that the hard over case ( i.e., one control 
surface locked at a certain deflection angle) shows greater excursion 
in the normal acceleration and roll angle with large values being 
achieved, than the other type of impairments. 

(3) It was showed that the impaired missile was able to track the 
unimpaired missile as long as there was another control surface to 
provide the appropriate corrective action and the control power of the 
failed surface. 

(4) A discussion was given of a possible method for transitioning from 
normal control laws to a set of reconfigured control lows for a 
missile with failed surface. The strategy requires the system to 
detect and isolate a surface failure and then transfer control to an 
off-line flight control computer which has been programmed with the 
new control laws (corrective actions). 

(5) The state-feedback autopilot was introduced as additional dynamic 
design in order to implement control to the unimpaired system. This 
design, analyzed in terms of its transient response, was found to 
meet the desired requirements. 
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(6) The performance of the coupled pitch and roll channel was found to be 
satisfactory and the unimpaired system proved to be robust. 



B. RECOMMENDATIONS 

Finally the following recommendations should be taken Into 
consideration: 

(1) More complex and other type of Impairments could be Investigated 
using the designed autopilot. In order to examine the overall 
controllability of the missile. 

(2) In order to make the system simpler, some of the returning gain loops 
of the state-feedback design must be eliminated. A further 
Investigation must then be conducted In order to examine If the 
performance of the resulting design remains robust. 

(3) Further analysis should be conducted If robustness were required In 
both the Input and output case simultaneously. 
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APPENDIX A 



PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
PITCH CHANNEL AUTOPILOT (UNIMPAIRED SYSTEM) 
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0 . 00000*00 


0 . 17890*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10000*01 


- 0 . 70190*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 66360*02 


- 0 . 16920*00 


- 0 . 10630*06 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

- 0 . 96200*02 


0 . 00000*00 

0 . 96200*02 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 18860*03 


- 0 . 96200*02 


0 . 96200*02 


0 . 23060*00 

0 . 00000*00 


0 . 11150*01 

0 . 00000*00 


0 . 17030*00 

0 . 00000*00 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


0 . 10090*02 


- 0 . 16300*00 


0 . 00000*00 


- 0 . 23060*00 

0 . 00030*00 


- 0 . 11150*01 

0 . 00000*00 


- 0 . 17030*00 

0 . 00000*00 


0 . 11150*01 


0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10090*02 


0 . 00000*00 


- 0 . 16300*00 


0 . 23060*00 

- 0 . 16300*00 


0 . 11150*01 

0 . 00000*00 


0 . 17030*00 

0 . 00000*00 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


0 . 10090*02 


0 . 00000*00 


0 . 00000*00 


- 0 . 23060*00 

0 . 00000*00 


- 0 . 11150*01 

- 0 . 16300*00 


- 0 . 17030*00 

0 . 00000*00 


0 . 11150*01 


0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10090*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

0 . 66320*03 


0 . 00000*00 

- 0 . 66320*03 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10000*01 


0 . 88710*03 


0 . 66320*03 


- 0 . 66320*03 



OPEN LOOP EIGENVALUES DET(SI-F)., 



1-5972AE+02, 1.89922E+01:-8.29043E+00. 8 



05911E+00S-3.75933E+00, 2.51AA1E+00 
65991E+00;-l. A3000E-01 :-l . A3000E-01i 



THE CONTROL DISTRIBUTION MATRIX G. 



0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 .9A20D+02 
0 . 1A30D+00 
O.OOOOD+00 
O.OOOOD+00 
0 . OOOOD+00 
0. AA32D+03 



O.OOOOD+00 
O.OOOOD+00 
0 . OOOOD+00 
0. OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
O.OOOOD+00 
0 . 9A20D+02 
0 . OOOOD+00 
-0.1A30D+00 
O.OOOOD+00 
O.OOOOD+00 
-0. AA32D+03 



O.OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
O.OOOOD+00 
O.OOOOD+00 
0. OOOOD+00 
0. OOOOD+00 
-0.9A20D+02 
0 . OOOOD+00 
0 . OOOOD+00 
-0.1A30D+00 
O.OOOOD+00 
0. AA32D+03 



O.OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0. OOOOD+00 
O.OOOOD+00 
0 . OOOOD+00 
O.OOOOD+00 
0.9A20D+02 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
-0 , 1A30D+00 
-0.AA32D+03 



0 . OOOOD+00 
-0.1353D+01 
0.6572D+01 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
0 . OOOOD+00 
-0 .230AD+00 
0 .230AD+00 
-0.230AD+00 
0.230AD+00 
0 . OOOOD+00 
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THE OPTIMAL FEEDBACK GAIN CONTROL MATRIX ... C = BINVXGTXS . . 



1 .A313D-01 
9.3121D-02 
-1.A313D-01 
-9.3121D-02 
1 . A313D-01 
9.3121D-02 
-1 .A313D-01 
-9.3121D-02 
-5.5193D-03 
-3.8930D-01 



-8.1227D-01 
-9.3121D-02 
8.1227D-01 
9 .3121D-02 
-8 .1227D-01 
-9.3121D-02 
8 .1227D-01 
9.3121D-02 
-2.3641D+00 
3.8930D-01 



-1.5A52D-01 
-8.17A1D-01 
1 .5A52D-01 
8.17A1D-01 
-1 .5A52D-01 
-8.17A1D-01 
1 .5A52D-01 
8 .17A1D-01 
-A.3398D-01 
9.2A0AD-02 



1 .570AD-02 
-1 .570AD-02 
1 .570AD-02 
-1 .570AD-02 
-2.A758D-01 



3.2362D-03 

-3.2362D-03 

3.2362D-03 

-3.2362D-03 

^.0967D-02 



-2. AA26D-01 
2.4A26D-01 
- 2 .^^ 4260-01 
2.4426D-01 
1.7860D-01 



C-LOOP OPTIMAL REG. E-VALUES ... DET( SI-F+G)CC) . . 



-9. 09132E+02:-2. 11193E+02:-l .78501E+01 : -9 . 79306 E+00 ; -5 . 22654E+00 : -6 . 87901 E-01 : 
-7 .46000E-02, 6 . 65991 E+00 : -1 . 46 233E-01 : -5 . 42351 E-06 : -1 . 43000 E-01 ; -1 . 43000 E-01 : 
-1 .43000E-01 : 
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APPENDIX B 



PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
ROLL CHANNEL AUTOPILOT (UNIMPAIRED SYSTEM) 





THE SYSTEM 


MATRIX "F"- 


•MATRIX . . . 






0.00000 

0.00000 


0.00000 

0.00000 


. 0.00000 
0.00000 


0.00000 

8.03A62 


0.00000 


0 . 00000 


1.00000 

0.00000 


0.00000 

0.00000 


0.00000 
0 . 00000 


0.00000 

0.00000 


0.00000 


0.00000 


0.00000 

0.00000 


-17.60000 

0.00000 


-8.00000 

0.00000 


0.00000 

0.00000 


0.00000 


0.00000 


-0.75325 

0.00000 


-0.88352 

0.00000 


0.35165 

0.00000 


-5.00000 
-0. A033A 


0.00000 


0.00000 


0.00000 

0.00000 


0.00000 

0.00000 


0.00000 

0.00000 


0.00000 

0.68291 


-6.00000 


0.00000 


-0.10271 

0.00000 


-0.120^8 

0.00000 


0.0A795 

0.00000 


lA. 31820 
-0.1A067 


-1A.1818A 


-15.00000 


0.10271 

-15.00000 


0.120^8 

0.00000 


-0 . 0A795 
0.00000 


-lA. 31820 
0.1A067 


1A.1818A 


0 . 00000 


0.10271 

0.00000 


0.120^8 

-15.00000 


-0.0A795 

0.00000 


-lA. 31820 
0.1A067 


1A.1818A 


0.00000 


-0.10271 

0.00000 


-0.120^48 

0.00000 


0.0A795 

-15.00000 


lA. 31820 
-0.1A067 


-1A.1818A 


0.00000 


0 . 00000 
-2698.83000 


0 . 00000 
-2698.83000 


0.00000 

2698.83000 


0 . 00000 
-188 . AOOOO 


0.00000 2698.83000 


OPEN LOOP 


EIGENVALUES. . 








.85057E+02:- 

.18517E+00:- 


9.27<462E+00, 
1 .50000E+01 


2.92679E+01 s 
■1 .50000E+01 : 


-2.35821E+00, 
-1 .50000E+01 : 


2.70296E+00 


:-8.89209E+00 
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THE CONTROL DISTRIBUTION MATRIX "G"-MATRIX . . 



0.00000 
0.00000 
0 . 00000 
0.00000 
0.00000 
-15.00000 
0 . 00000 
0.00000 
0.00000 
2698.83000 



0 . 00000 
0 . 00000 
0.00000 
0.00000 
0 . 00000 
0 . 00000 
-15.00000 
0 . 00000 
0 . 00000 
-2698.85000 



0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

-15.00000 

0.00000 

-2698.83000 



0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0 . 00000 
0.00000 
-15 . 00000 
2698.83000 



0.00000 
0.00000 
17.60000 
0.88352 
0 . 00000 
0.120A8 
-0.120A8 
-0.12048 
0.12048 
0 . 00000 



THE OPTIMAL FEEDBACK GAIN CONTROL MATRIX ... C=BINVXGTX$ .. . 



5.0881D-01 

-5.0881D-01 

-5.0881D-01 

5.0881D-01 

1.7841D-01 



5.6129D-01 

-5.6129D-01 

-5.6129D-01 

5.6129D-01 

-1.5987D-01 



-3.0819D-02 
3.0819D-02 
3. 0819D-02 
-5.0819D-02 
1.1143D-01 



-1.1772D+00 

1.1772D+00 

1.1772D+00 

-1.1772D+00 

3.8252D+00 



1.0775D+00 
-1.0775D+00 
-1 . 0775D+00 
1 . 0775D+00 
-3.4606D+00 



-2.7569D-01 
-1 .3852D-01 
-1 .3852D-01 
1 .3852D-01 
5.6129D-01 



C-LOOP OPTIMAL REG. E-VALUES ... D£T(SI-F+GXC) . . 

-5.40070E+03:-5.99206E+00, 5. 34991 E+ 00: -9. 29732E+00 s-6 .57672E+00, 8 . 26712E-01 
-1 .45396 E+00: -2. 12132E+01 : -2.121 32E+01 : -2.1 2132E+01 : 
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APPENDIX C 



PLANT SYSTEM AND INPUT MATRICES OF CONTINUOUS OPEN LOOP 
COUPLED PITCH AND ROLL CHANNELS AUTOPILOTS 
(UNIMPAIRED SYSTEM) 



THC CONTROL DISTRIBUTION MATRIX 0.. 



0.0000 D «00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


O . OOOOD^OO 


0.00000*00 


0.00000*00 


0.00000*00 


-0.13530*01 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


O . OOOOD^OO 


0.00000*00 


0.00000*00 


0.00000*00 


0.65720*01 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.0000 D «00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000^00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.0000 D «00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


O . OOOOD^OO 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


■0.9420 D ^02 


0.94200*02 


-0.94200*02 


0.94200*02 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


■0.1430 D «00 


0.00000*00 


0.00000*00 


0.00000*00 


-0.23040*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


O . OOOOD^OO 


-0.14300*00 


0.00000*00 


0.00000*00 


0.23040*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0. 14300*00 


0.00000*00 


-0.23040*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0.14300*00 


0.23040*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.44320*03 


-0.44320*03 


0.44320*03 


-0.44320*03 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.17600*02 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.88350*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0.15000*02 


0.00000*00 


0.00000*00 


0.00000*00 


0.12050*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0. 15000*02 


0.00000*00 


0.00000*00 


-0. 12050*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0.15000*02 


0.00000*00 


-0.12050*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


-0.15000*02 


0.12050*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.00000*00 


0.26990*04 


-0.26990*04 


-0.26990*04 


0.26990*04 


0.00000*00 



OPEN LOOP CIGEMVALUeS DET(SI-F).. 



-1 .59724 E *02, 1 . 89922 E *0 1 x -8 .2904 JE *0 0 . 8 .0591 IE *00 x -3 . 75933 E »00 , 2.51441 E »00 x 
- 1 .43935 E -01 :-5 .3141 7 E -06 x - 7 . 46000 E -02 t 6 . 65 99 IE * 00 1 .43000 E -01 1 .43000 E -0 1 x 
-1.43000 E -01 :- l .85057 E *02:-9.27462 E *00 . 2.92679 E + 01 x -2 . 3582 1 E * 00 . 2.70296 E ^00 x 
“8.89209 E *00 :-5 . 185 1 7 E ^00 x - 1 . 50000 E *0 1 x - 1 .50000 E ^01 1 .50000 E *01 x 
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OPEN LOOP dynamics matrix 



F 



• 0 . 15000*03 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


o . ooooo*oo 


0 . 00000*00 


0 . 15000*03 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 13530*01 


0 . 00000*00 


0 . 10000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















■ 0 . 65720*01 


0 . 00000*00 


- 0 . 50000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10000*01 


0 . 00000*00 


0 . 00000*00 


- 0 . 59260*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 44360*02 


- 0 . 14920*00 


0 . 00000*00 


0 . 00000*00 


0 . 17890*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10000*01 


- 0 . 70190*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 44360*02 


- 0 . 14920*00 


- 0 . 10430*04 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00001^*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 18840*03 


- 0 . 94200*02 


0 . 94200*02 


■ 0 . 9 < 20D*02 


0 . 94200*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 230AD*00 


0 . 11150*01 


0 . 17030*00 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


0 . 10090*02 


- 0 . 14300*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















• 0 . 23040*00 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 11150*01 


0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10090*02 


0 . 00000*00 


- 0 . 14300*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 23040*00 


0 . 11150*01 


0 . 1703D*Oa 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


0 . 10090*02 


0 . 00000*00 


0 . 00000*00 


■ 0 . 14300*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















• 0 . 23040*00 


- 0 . 11150*01 


- 0 . 17030*00 


0 . 11150*01 


0 . 17030*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10090*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 14300*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


O . O000D * 0O 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10000*01 


0 . 88710*03 


0 . 44320*03 


- 0 . 44320*03 


0 . 44320*03 


- 0 . 44320*03 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 : 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 80350*01 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 17600*02 


- 0 . 80000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 75330*00 


- 0 . 88350*00 


0 . 35170*00 


- 0 . 50000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 40330*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 60000*01 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 68290*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 10270*00 


- 0 . 12050*00 


0 . 47950-01 


0 . 14320*02 


- 0 . 14180*02 


- 0 . 15000*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


- 0 . 14070*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10270*00 


0 . 12050*00 


- 0 . 47950-01 


- 0 . 14320*02 


0 . 14180*02 


0 . 00000*00 


- 0 . 15000*02 


0 . 00000*00 


0 . 00000*00 


0 . 14070*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 10270*00 


0 . 12050*00 


- 0 . 47950-01 


- 0 . 14320*02 


0 . 14180*02 


0 . 00000*00 


0 . 00000*00 


• 0 . 15000*02 


0 . 00000*00 


0 . 14070*00 
















0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

- 0 . 15000*02 


0 . 00000*00 

- 0 . 14070*00 


- 0 . 10270*00 


- 0 . 12050*00 


0 . 47950-01 


0 . 14320*02 


- 0 . 14180*02 


0 . 00000*00 


0 . 00000*00 


0 . 00000*00 

0 . 00000*00 

0 . 26990*04 


0 . 00000*00 

0 . 00000*00 

0 . 26990*04 


0 . 00000*00 

0 . 00000*00 

- 0 . 18840*03 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 00000*00 


0 . 00000*00 

0 . 26990*04 


0 . 00000*00 

- 0 . 26990*04 
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APPENDIX D 



COMPUTER PROGRAMS 

In order to plot the state variables of interest and compare them, the 
following computer programs were developed in WATFIV and used, in order 
to change the format of the OPTPLOT DATA, into a fixed format data file, 80 
characters in length with space between the values, to be used as input file 
for the EASYPLOT program. 

1. Uncoupled Pitch Channel and for q . Sp and ty 



1 


$JOB 


REAL A,B,C,D 


2 




N=1 


3 




WHILE (N.LE.lOO) DO 


4 




READ(5,10)A 


5 




READ(5,15)B 


6 




READ(5,20)C 


7 




READ(5,25)D 


8 




WRITE(6,30)A,B,C,D 


9 


10 


F0RMAT(E14.7) 


10 


15 


F0RMAT(56X,E14.7) 


11 


20 


F0RMAT(A2X, E14.7) 


12 


25 


F0RMAT(42X,E1A.7) 


13 


30 


F0RMAT(8F15.8) 


14 




N = N+1 


15 




END WHILE 


16 




STOP 


17 




END 



GENTRY 



2. Uncoupled Roll Channel end for P. and Sp 



1 


$JOB 


REAL A,B,C,D,E 


2 




N = 1 


3 




WHILE (N.LE.lOO) DO 


4 




READ(5,10)A 


5 




READ(5,15)B,C 


6 




READ(5,20)D 


7 




READ(5,25)E 


8 




WRITE(6,50)A,B,C,E 


9 


10 


FORMAT (El < 4 . 7) 


10 


15 


F0RMAT(14X,2E1A.7) 


11 


20 


FORMAT (El A. 7) 


12 


25 


F0RMAT(E1A.7) 


13 


30 


F0RMAT(4F15.8) 


14 




N = N+1 


15 




END WHILE 


16 




STOP 


17 




END 



GENTRY 
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3. Coupled Pitch and Roll Channels and fQr_g . &> andJb 



$JOB 



1 




REAL A,B,C,D,E,L,G 


2 




N = 1 


3 




HHILE (N.LE.lOO) DO 






READ(5',10)A 


5 




READ(5,15)B 


6 




READ(5,20)C 


7 




READ(5,25)D 


8 




READ(5,30)E 


9 




READ(5,35)L 


10 




READ(5,A0)G 


11 




WRITE(6,<i5)A,C,D,E 


12 


10 


F0RMAT(E1A.7) 


13 


15 


F0RMAT(E1A.7) 


1^ 


20 


F0RMAT(56X,E1A.7) 


15 


25 


F0RMAT(A2X,E1A.7) 


16 


30 


F0RMAT(<«2X,El<i .7) 


17 


35 


F0RMAT(E1A.7) 


18 


^0 


F0RMAT(E1A.7) 


19 


^5 


F0RMAT(AF15.8) 


20 




N = N+1 


21 




END WHILE 


22 




STOP 


23 




END 



$ENTRY 



4. Coupled Pitch and Roll Channels end for P. ■!> and &> 



$JOB 



1 




REAL A,B,C,D,E,F,G 


2 




N = 1 


3 




WHILE (N.LE.lOO) DO 






READ(5,10)A 


5 




READ(5,15)B 


6 




READ(5,20)C 


7 




READ(5,25)D 


8 




READ(5,30)E 


9 




READ(5,35)F 


10 




READ(5,40)G 


11 




WRITE(6,45)A,E,F,G 


12 


10 


F0RMAT(E1<4.7) 


13 


15 


F0RMATCE14.7) 


1^ 


20 


F0RMATCE1A.7) 


15 


25 


F0RMAT(E1A.7) 


16 


30 


F0RMAT(56X,E1A.7) 


17 


35 


F0RMAT(E1A.7) 


18 


^0 


F0RMAT(42X,E1A.7) 


19 


^5 


F0RMAT(4F15.8) 


20 




N=N+1 


21 




END WHILE 


22 




STOP 


25 




END 



GENTRY 
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APPENDIX E 



COUPLED STATE-FEEDBACK DESIGN INPUT DATA 
FOR POPLAR PROGRAM 
(UNIMPAIRED SYSTEM) 



1 1 

100.0 0.1 0.1 
1.0 0.00 0.0 

0.5 0.6 .1 0 

0 10 0 5 S 1000 0 

232325102325102S2S 

-0. 1500000000000D*05 O.OOOOOOOOOOOOOD^OO 
0 . 0000000000000 D «00 o.oooooooooooooD«oo 
O.OOOOOOOOOOOOOD^OO 0.00000000000000*00 
0. 15OOOOOOOOOO0D«OS 0. OOOOOOOOOOOOOD^OO 
O.OOOOOOOOOOOOOD^OO 0.000000000 00 OOD^OO 



O.OOOOOOOOOOOOOD^OO 

O.OOOOOOOOOOOOOD^OO 

o.oooooooooooooD«oo 

0.0000000000000D«00 
O.OOOOOOOOOOOOOD^OO 
0 . 00000000000000*00 
0 . 00000000000000*00 
-0.50000000000000*01 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
-0.59260000000000*02 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0.00000000000000*00 
-0. U 920000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
-0.16920000000000*00- 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.0000000000000D*00 
-0.18860000000000*03- 
0.96200000000000*02 
0 . 00000000000000*00 
0.00000000000000*00 
0.23062000000000*00 
-0. 1703000000000O*00 
-0. 16300000000000*00 



0 . 00000000000000*00 
0 . 10000000000000*01 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000 * 00 * 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooo«ooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooD*oo 
0 . 10000000000000*01 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooD*oo 
0.00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 10000000000000 * 01 - 
0.00000000000000*00 
0.00000000000000*00 
0. 00000000000000*00 
0.00000000000000*00 
o.oooooooooooooo*oo 

-0. l063S60000000O*06 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0 . 00000000000000*00 
0.00000000000000*00 
0.96200000000000*02 
0.00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.11156000000000*01 
0 . 00000000000000*00 
0.0000000000000O*00 



0 . 00000000000000*00 
o.ooooooooooo^oo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0.00000000000000*00 
0 . 00000000000000*00 
-0.65720000000000*01 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000 * 00 - 
0.00000000000000*00 
o.oooooooooooooo*oo 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
>0.70190000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0. 0000000000000O*00- 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.96200000000000*02- 
0.00000000000000*00 
0.00000000000000*00 
0.00000000000000*00 
0.17030000000000*00- 
0.00000000000000*00 
0.00000000000000*00 



0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.13531000000000*01 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 
o.oooooooooooooo*oo 

0.0000000000000D*00 

0 . 00000000000000*00 

0.00000000000000*00 

o.oooooooooooooD*oo 

>0.66360000000000*02 

0.17887000000000*01 

0.00000000000000*00 

0.00000000000000*00 

0.00000000000000*00 

0.00000000000000*00 

o.oooooooooooooD*oo 

0.00000000000000*00 

0.00000000000000*00 

0.00000000000000*00 

o.oooooooooooooD*oo 

0 . 00000000000000*00 

0.66360000000000*02 

0.0000000000000O*00 

0.00000000000000*00 

0.00000000000000*00 

o.oooooooooooooo*oo 

0.00000000000000*00 

0.00000000000000*00 

0.9620000000000D*02 

0.00000000000000*00 

0 . 00000000000000*00 

o.oooooooooooooo*oo 

0.1 1 1S6000000000*01 
0.10092000000000*02 
0 . 00000000000000*00 
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00000000000000*00 
00000000000000* 00 
00000000000000*00 
iii 5 <;ooooooooo*oi' 
00000000000000*00 
14300000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
17 O 3 OOOO 00000 O * 0 O - 
00000000000000*00 
14300000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000 * 00 - 
11154000000000*01 
10092000000000*02 
14300000000000*00 
. 00000000000000*00 
. 00000000000000*00 
. 00000000000000*00 
. 00000000000000*00 
44321100000000 * 03 - 
00000000000000*00 
. 00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
oooooooooooooo*oo 
oooooooooooooo*oo' 
00000000000000*00 
8 OOO 000000000 O *01 
00000000000000*00 
00000000000000*00 
O 0000000000000*00 
00000000000000*00 
00000000000000 * 00 - 
5000000000000 D *01 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
60000000000000*01 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
12048000000000*00 
15 OOOOOO 0000 OO * O 2 
14067000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
47950000000000 - 01 - 
15000000000000*02 
00000000000000*00 
00000000000000*00 



0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
- 0 . 17030000000000*00 0 . 

o.oooooooooooooo*oo-o. 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
- 0 . 11154000000000 * 01 - 0 . 
0 . 10092000000000*02 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
- 0 . 23042000000000 * 00 - 0 . 
0 . 17030000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
o.oooooooooooooo*og-o. 

- 0 . 44321100000000*03 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 80346200000000*01 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
O . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
O . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
■ 0 . 75325000000000 * 00-0 
0 . 00000000000000*00 0 . 
0 . 00000000000000 * 00 - 0 . 
0 . OO 00 O 00000000 O *00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 68291000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 47950000000000-01 0 . 
O . OOOOOOOO 000 OOO *00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 14318200000000*02 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 
0 . 00000000000000*00 0 . 



00000000000000*00 0 
00000000000000*00 0 
oooooooooooooo*oo-o 

11154000000000*01 0 
10092000000000*02 0 
00000000000000*00 0 
00000000000000*00 0 
00000000000000*00 0 
23042000000000*00 0 
17030000000000*00 0 
00000000000000*00 0 
0 OO 0 OOOOO 0000 O *00 0 
00000000000000*00 0 
00000000000000*00 0 
11154000000000 * 01-0 
00000000000000*00 0 
00000000000000*00 0 
. 00000000000000*00 0 
. 00000000000000*00 0 
. 00000000000000*00 0 
. 00000000000000*00 0 
. 10000000000000*01 0 
, 44321100000000 * 03-0 
. 00000000000000*00 0 
. 00000000000000*00 0 
00000000000000*00 0 
. 00000000000000*00 0 
00000000000000*00 0 
00000000000000*00 0 
. 00000000000000*00 0 . 
oooooooooooooo*oo 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
oooooooooooooo*oo o. 
00000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
O 0000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000 * 00 - 0 . 
O 0000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
OOO 000 OO 00000 O *00 0 . 
00000000000000*00 0 . 
88352000000000*00 0 . 
00000000000000*00 0 . 
40334000000000*00 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
oooooooooooooo*oo 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
oooooooooooooo*oo 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000 * 00 - 0 . 
14318200000000 * 02 - 0 . 
O 0000000000000*00 0 . 
00000000000000*00 0 . 
OOOOOOOOOOOOOO * O 0 0 . 
00000000000000*00 0 . 
102710 OO 00000 O *00 0 . 
14181840000000*02 0 . 
00000000000000*00 0 . 
00000000000000*00 0 . 
00000000000000*00 0 



00000000000000*00 
00000000000000*00 
23042000000000*00 
17030000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
11154000000000*01 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
17030000000000*00 
00000000000000*00 
00000000000000*00 
. OOO 0000 O 00000 O *00 
00000000000000*00 
. 00000000000000*00 
. 00000000000000*00 
. 88712400000000*03 
44321100000000*03 
. 00000000000000*00 
. 00000000000000*00 
00000000000000*00 
. 00000000000000*00 
. 00000000000000*00 
. 00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
10000000000000*01 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
17600000000000*02 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
35165000000000*00 
00000000000000*00 
00000000000000*00 
oooooooooooooo*oo 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
00000000000000*00 
oooooooooooooo*oo 

lO 271 OOOOOO 0 OO *00 
14181840000000*02 
00000000000000*00 
00000000000000*00 
oooooooooooo0O*00 
00000000000000* 00 
12048000000000*00 
00000000000000*00 
14067000000000*00 
00000000000000*00 
00000000000000*00 
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0 . 00000000000000*00 
0 . 00000000000000*00 
-0.14318200000000*02 
-0. lS00000000000O*02 
0 . 00000000000000*00 
0 . 00000000000000*00 
o.oooooooooooooo*oo 

-0.10271000000000*00- 
-0. 14181840000000*02 
-0. 15000000 OOOOOO* 02- 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0 . 00000000000000*00 
0.26988300000000*04- 
-0.18840000000000*03 

o.oooooooooooooo*oo 

0 . 00000000000000*00 

0.00000000000000*OQ 

o.oooooooooooooo*oo 

0 . 00000000000000*00 

0 . 00000000000000*00 

0.65720000000000*01 

0 . 00000000000000*00 
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